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Abstract 

These lectures provide a modern introduction to selected topics in the physics 
of ultrarelativistic heavy ion collisions which shed light on the fundamental 
theory of strong interactions, the Quantum Chromodynamics. The emphasis is 
on the partonic forms of QCD matter which exist in the early and intermediate 
stages of a collision — the colour glass condensate, the glasma, and the quark- 
gluon plasma — and on the effective theories that are used for their description. 
These theories provide qualitative and even quantitative insight into a wealth 
of remarkable phenomena observed in nucleus-nucleus or deuteron-nucleus 
collisions at RHIC and/or the LHC, like the suppression of particle production 
and of azimuthal correlations at forward rapidities, the energy and centrality 
dependence of the multiplicities, the ridge effect, the limiting fragmentation, 
the jet quenching, or the dijet asymmetry. 



*Based on lectures presented at the 2011 European School of High-Energy Physics (ESHEP2011), 7-20 September 2011, 
Cheile Gradistei, Romania. 
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1 Introduction 

With the advent of the high-energy colliders RHIC (the Relativistic Heavy Ion Collider operating at 
RHIC since 2000) and the LHC (the Large Hadron Collider which started operating at CERN in 2008), 
the physics of relativistic heavy ion collisions has entered a new era: the energies available for the 
collisions are high enough — up to 200 GeV per interacting nucleon pair at RHIC and potentially up 
to 5.5 TeV at the LHC (although so far one has reached 'only' 2.76 TeV) — to ensure that new forms 
of QCD matter, characterized by high parton densities, are being explored by the collisions. These new 
forms of matter refer to both the wavefunctions of the incoming nuclei, prior to the collision, which 
develop high gluon densities leading to colour glass condensates, and the partonic matter produced in 
the intermediate stages of the collision, which is expected to form a quark-gluon plasma. The asymptotic 
freedom property of QCD implies that these high-density forms of matter are weakly coupled (at least 
in so far as their bulk properties are concerned) and hence can be studied via controlled calculations 
within perturbative QCD. But such studies remain difficult and pose many challenges to the theorists: 
precisely because of their high density, these new forms of matter are the realm of collective, non- 
linear phenomena, whose mathematical description often transcends the ordinary perturbation theory. 
Moreover, there are also phenomena (first revealed by the experiments at RHIC) which seem to elude a 
weak-coupling description and call for non-perturbative techniques. 

These challenges stimulated new ideas and the development of new theoretical tools aiming at a 
fundamental understanding of QCD matter under extreme conditions : high energy, high parton densi- 
ties, high temperature. The ongoing experimental programs at RHIC and the LHC provide a unique and 
timely opportunity to test such new ideas, constrain or reject models, and orient the theoretical develop- 
ments. Over the last decade, the experimental and theoretical efforts have gone hand in hand, leading 
to a continuously improving physical picture, which is by now well rooted in QCD. The purpose of 
these lectures is to provide an introduction to this physical picture, with emphasis on those aspects of 
the dynamics for which we are confident to have a reasonably good (although still far from perfect) un- 
derstanding from first principles, i.e. from the Lagrangian of Quantum Chromodynamics. These aspects 
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Fig. 1: Schematic representation of the various stages of a HIC as a function of time t and the longitudinal 
coordinate z (the collision axis). The 'time' variable which is used in the discussion in the text is the proper time 
r = \/t 2 — z 2 , which has a Lorentz-invariant meaning and is constant along the hyperbolic curves separating 
various stages in this figure. 



concern the partonic stages of a heavy ion collision, at sufficiently early times. These are also the stages 
to which refers most of the experimental and theoretical progress over the last decade. 

2 Stages of a heavy ion collision: the case for effective theories 

The theoretically motivated space-time picture of a heavy ion collision (HIC) is depicted in Fig.[T] This 
illustrates the various forms of QCD matter intervening during the successive phases of the collision: 

1. Prior to the collision, and in the center-of-mass frame (which at RHIC and the LHC is the same as 
the laboratory frame), the two incoming nuclei look as two Lorentz-contracted 'pancakes', with a 
longitudinal extent smaller by a factor 7 ^ 100 (the Lorentz boost factor) than the radial extent in 
the transverse plane. As we shall see, these 'pancakes' are mostly composed with gluons which 
carry only tiny fractions x <C 1 of the longitudinal momenta of their parent nucleons, but whose 
density is rapidly increasing with 1/x. By the uncertainty principle, the gluons which make up 
such a high-density system carry relatively large transverse momenta. A typical value for such a 
gluon in a Pb or Au nucleus is k± ~ 2 GeV for x = 10 -4 . By the 'asymptotic freedom' property of 
QCD, the gauge coupling which governs the mutual interactions of these gluons is relatively weak. 
This gluonic form of matter, which is dense and weakly coupled, and dominates the wavefunction 
of any hadron (nucleon or nucleus) at sufficiently high energy, is universal — its properties are the 
same form all hadrons. It is known as the colour glass condensate (CGC). 

2. At time r = 0, the two nuclei hit with each other and the interactions start developing. The 
'hard' processes, i.e. those involving relatively large transferred momenta Q > 10 GeV, are those 
which occur faster (within a time r ~ l/Q, by the uncertainty principle^]). These processes are 
responsible for the production of 'hard particles', i.e. particles carrying transverse energies and 
momenta of the order of Q. Such particles, like (hadronic) jets, direct photons, dilepton pairs, 
heavy quarks, or vector bosons, are generally the most striking ingredients of the final state and 
are often used to characterize the topology of the latter — e.g., one speaks about 'a dijet event', cf. 
Fig.[2]left, or 'a photon-jet' event, cf. Fig.[2]right. 

3. At a time r ^ 0.2 fm/c, corresponding to a 'semi-hard' transverse momentum scale Q ~ 1 GeV, 

throughout these notes, we shall generally use the natural system of units h — c — ks — ^, so in particular there is no 
explicit factor h in the uncertainty principle. Yet, in some cases, we shall restore this factor for more clarity. 
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the bulk of the partonic constituents of the colliding nuclei (meaning the gluons composing the 
respective CGCs) are liberated by the collision. This is when most of the 'multiplicity' in the 
final state is generated; that is, most of the hadrons eventually seen in the detectors are produced 
via the fragmentation and the hadronisation of the initial-state gluons liberated at this stage. But 
before ending up in the detectors, these partons undergo a complex evolution. Just after being 
liberated, they form a relatively dense medium, whose average density energy in Pb+Pb collisions 
at the LHC is estimated as e > 15 GeV/fm 3 ; this is about 10 times larger than the density of 
nuclear matter and 3 times larger than in Au+Au collisions at RHIC. This non-equilibrium state of 
partonic matter, which besides its high density has also other distinguished features to be discussed 
later, is known as the glasma. 
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Fig. 2: A couple of di-jets events in Pb+Pb collisions at ATLAS (left) and CMS (right). 



4. If the produced partons did not interact with each other, or if these interactions were negligible, 
then they would rapidly separate from each other and independently evolve (via fragmentation 
and hadronization) towards the final-state hadrons. This is, roughly speaking, the situation in 
proton-proton collisions. But the data for heavy ion collisions at both RHIC and the LHC ex- 
hibit collective phenomena (like the 'elliptic flow' to be discussed later) which clearly show that 
the partons liberated by the collision do actually interact with each other, and quite strongly. A 
striking consequence of these interactions is the fact that this partonic matter rapidly approaches 
towards thermal equilibrium : the data are consistent with a relatively short thermalization time, 
of order r ~ 1 fm/c. This is striking since it requires rather strong interactions among the partons, 
which can compete with the medium expansion: these interactions have to redistribute energy and 
momentum among the partons, in spite of the fact that the latter separate quite fast away from 
each other. Such a rapid thermalization seems incompatible with perturbative calculations at weak 
coupling and represents a main argument in favour of a new paradigm: the dense partonic matter 
produced in the intermediate stages of a HIC may actually be a strongly coupled fluid. 

5. The outcome of this thermalization process is the high-temperature phase of QCD known as the 
quark-gluon plasma. The abundant production and detailed study of this phase is the Holy Grail 
of the heavy ion programs at RHIC and the LHC. The existence of this phase is well established 
via theoretical calculations on the lattice, but its experimental production within a HIC is at best 
ephemeral: the partonic matter keeps expanding and cooling down (which in particular implies 
that the temperature is space and time dependent, i.e. thermal equilibrium is reached only locally) 
and it eventually hadronizes — the 'coloured' quark and gluons get trapped within colourless 
hadrons. Hadronization occurs when the (local) temperature becomes of the order of the critical 
temperature T c for deconfinement, known from lattice QCD studies as T c ~ 150-1-180 MeV. In 
Pb+Pb collisions at the LHC, this is estimated to happen around a time r ^ 10 fm/c. 
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6. For larger times 10 < r < 20 fm/c, this hadronic system is still relatively dense, so it preserves lo- 
cal thermal equilibrium while expanding. One then speaks of a hot hadron gas, whose temperature 
and density are however decreasing with time. 

7. Around a time r ~ 20 fm/c, the density becomes so low that the hadrons stop interacting with 
each other. That is, the collision rate becomes smaller than the expansion rate. This transition 
between a fluid state (where the hadrons undergo many collisions) and a system of free particles is 
referred to as the freeze-out. From that moment on, the hadrons undergo free streaming until they 
reach the detector. One generally expects that the momentum distribution of the outgoing particles 
is essentially the same as their thermal distribution within the fluid, towards the late stages of 
the expansion, just before the freeze-out. This assumption appears to be confirmed by the data: 
the particle spectra as measured by the detectors can be well described as thermal (Maxwell- 
Boltzmann) distributions, with only few free parameters, like the fluid temperature and velocity at 
the time of freeze-out. This is generally seen as an additional argument in favour of thermalization, 
but one must be cautious on that, since the mechanism of hadronisation itself can lead to spectra 
which are apparently thermal. As a matter of fact, the freeze-out temperature extracted from the 
ratios of particle abundances at RHIC appears to be the same, Tf ~ 170 MeV, in both Au+Au and 
p+p collisions, while of course no QGP phase is expected in p+p. 

Although extremely schematic, this simple enumeration of the various stages of a HIC already 
illustrates the variety and complexity of the forms of matter traversed by the QCD matter liberated by 
the collision on its way to the detectors. In principle, all these forms of matter and their mutual transfor- 
mations admit an unambiguous theoretical description in the framework of Quantum Chromodynamics, 
which is the fundamental theory of strong interactions. But although this theory exists with us for about 
40 years, it is still far from having delivered all its secrets. Indeed, in spite of the apparent simplicity of its 
Lagrangian, which looks hardly more complicated than that of the quantum electrodynamics (the theory 
of photons and electrons), the QCD dynamics is considerably richer and more complicated — which is 
why it can accommodate so many phases! What renders the theoretical study of HIC's so difficult is the 
extreme complexity of the relevant forms of hadronic matter, characterized by high (parton or hadron) 
densities and strong collective phenomena. For a theorist, the most efficient way to try and organize this 
complexity is to build effective theories. 

An 'effective theory' should not be confused with a 'model': its main purpose is not to provide 
a heuristic description of the data using some physical guidance together with a set of free parameters. 
Rather, it aims at a fundamental understanding and its construction is always guided by the underlying 
fundamental theory — here, QCD. Specifically, an effective theory is a simplified version of the fun- 
damental theory which includes the 'soft' (i.e. low energy and momentum) degrees of freedom (d.o.f.) 
required for the description of the physical phenomena occurring at a relatively large space-time scale, 
but ignores the 'hard' d.o.f. with higher energies and momenta. More precisely, the hard modes cannot 
be totally ignored — they interact with the soft modes and thus affect the properties of the effective the- 
ory — , rather they are 'integrated out' via some coarse-graining (or 'renormalization group') procedure, 
which can be perturbative or non-perturbative. 

If the coupling is weak (g <C 1), the 'hard-soft' interactions can be treated in perturbation theory 
and then the effective theory emerges as a controlled approximation to the original theory. This gen- 
erally amounts to computing Feynman graphs with hard loop momenta and soft external legs. By the 
uncertainty principle, the hard modes are localized on short space-time distances, so their net effect is to 
provide quasi-local vertices, or 'parameters' — like effective masses and couplings — in the effective 
Lagrangian for the soft modes. But even at weak coupling, one often has to deal with a large, or even 
infinite, number of Feynman graphs at any given order in g, because the contributions due to individual 
graphs are enhanced by the large disparity of scales between the hard and soft d.o.f. and/or by the high 
density of medium constituents. This is where the effective theory is most useful: it allows us to 're- 
sum' (modulo some approximations) a large number of Feynman graphs of the original field theory and 
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replace their effects by a small number of 'parameters' in the effective Lagrangian. 

When the coupling is relatively strong, g > 1, standard perturbation theory (the expansion in 
powers of g) is bound to fail and the construction of effective theories becomes more problematic. So 
long as the coupling is just moderately strong, say g ~ 0(1), there is still hope that some insightful 
resummations of the perturbation theory, as based on the proper identification of the relevant d.o.f., 
may reasonably work — we shall later encounter some examples in that sense. If, in some regime, the 
coupling happens to be even stronger, perturbation theory brings no guidance anymore, and there is no 
systematic method to construct effective theories. They can merely be postulated on the basis of general 
physical considerations, like the symmetries of the fundamental theory. In such a case, the effective 
masses or coupling constants are generally treated as free parameters, to be matched against the data or, 
in some cases, against lattice QCD calculations. Effective theories may also emerge for rather deep and 
unexpected reasons, as we shall see on the example of the gauge/string duality later on. 

If in the previous discussion we mentioned both weak and strong coupling scenarios, is because 
in QCD — and indeed in any of the fundamental field theories in Nature — the coupling 'constant' is 
not fixed: it 'runs' with the typical momenta exchanged in the interactions, meaning that it is different 
when probing the physics on different space-time scales. What is essential about QCD is the property of 
asymptotic freedom : the fact that the coupling becomes weaker on shorter distances, or with increasing 
momentum transfer. Given our experience with electromagnetism, this property may look counterintu- 
itive. In QED, the electric charge of the nucleus inside an atom is well known to be screened by the 
surrounding electron cloud, so that the atom appears electrically neutral from far away. Similarly, the 
electric charge of an electron is screened by electron-positron virtual pairs which pump up from the 
vacuum, with the result that the effective charge « e . m .(^) = e 2 (i?)/47r decreases with the distance R 
from the electron. But in QCD, there is ^/-screening: the effective colour charge of a quark or gluon, 
as measured by its coupling 'constant' a s (R) = g 2 {R)/^ii, increases with R or with decreasing the 
transferred momentum Q. (Recall that Q ~ 1/R by the uncertainty principle.) Specifically, one has 
a s (Q 2 ) ~ 0.1 when Q = 100 GeV (the typical scale for electroweak physics and also for hadronic jets 
at the LHC). This 'asymptotic freedom' is, of course, the ultimate reason behind the success story of 
perturbative QCD in relation with 'hard' processes. It also justifies the use of perturbation theory for 
integrating out the 'hard' d.o.f. in the construction of effective theories for the 'soft' ones. 

But there is also the reverse of the medal: with decreasing Q below 100 GeV, the QCD coupling 
is increasing, albeit slowly, according to 



so that e.g. a s (Q 2 ) ~ 0.4 when Q = 2 GeV. Formally, Eq. ([!]) predicts that the coupling diverges when 
Q = Aqcd, but this equation cannot be trusted for Q < 1 GeV, as it has been obtained in perturbation 
theory. The fate of the QCD coupling for Q ~ Aqcd is still under debate, but various non-perturbative 
approaches suggest that a s (Q 2 ) should (roughly) saturate at a value close to one. For all purposes, this 
is very strong coupling (e.g. it corresponds to g ~ 3). 

After this digression through the general scope of an effective theory and the QCD running cou- 
pling, let us return to the main stream of our presentation, namely, the phases of QCD as probed in a HIC. 
Some key ideas, that will be succinctly mentioned here and developed in more detail in the remaining 
part of these lectures, are as follows: 

(i) The different stages of a HIC involve different forms of hadronic matter with specific active 
degrees of freedom. Their theoretical description requires different effective theories. 

(ii) During the early stages of the collision — the colour glass condensate and the glasma — the 
parton density is very high, the typical transverse momenta are semi-hard (a few GeV), and the QCD 
coupling is moderately weak, say a s ~ 0.3. In this case, perturbation theory is (at least, marginally) 
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valid, but it goes beyond a straightforward expansion in powers of a s . The construct the corresponding 
effective theory, one needs to resum an infinite class of Feynman graphs which are enhanced by high- 
energy and high gluon density effects. This has been done in the recent years, led to a formalism — 
the CGC effective theory — which offers a unified description from first principles for both the nuclear 
wavefunctions prior to the collision and the very early stages of the collision. A key ingredient in this 
construction is the proper recognition of the relevant d.o.f. : quasi-classical colour fields. The concept 
of field is indeed more useful in this high-density environnement than that of particle, since the phase- 
space occupation numbers are large 1), meaning that the would-be 'particles' overlap with each 
other and thus form coherent states, which are more properly described as classical field configurations. 

(iii) At later stages, the partonic matter expands, the phase-space occupation numbers decrease, 
and the concept of particle becomes again meaningful: the classical fields break down into particles. If 
these particles are weakly coupled (as one may expect by continuity with the previous stages), then their 
subsequent evolution can be described by kinetic theory. This is an effective theory which emerges under 
the assumption that the mean free path between two successive collisions is much longer than any other 
microscopic scale (like the duration of a collision or the Compton wavelength A = 1/ k± of a parti- 
cle). Over the last years, kinetic theory has been extensively derived from QCD at weak coupling, but 
the results appear to be deceiving: for instance, they cannot explain the rapid thermalization suggested 
by the data at RHIC and the LHC. (The thermalization times predicted by perturbative QCD are much 
larger, r > 10 fm/c.) Several alternative solutions have been proposed so far, but the final outcome is 
still unclear. One of these proposals is that the softer modes, which keep large occupation numbers and 
should be better described as classical fields, become unstable due to the anisotropy in the momentum 
distribution of the harder particles (which in turns follows from the disparity between longitudinal and 
transverse expansions). But numerical simulations of the coupled system soft fields-hard particles leads 
to thermalization times which are still too large. Another suggestion is that the partonic matter is (mod- 
erately) strongly coupled — the QCD coupling could indeed become larger, because of the system being 
more dilute. In such a scenario, a candidate for an effective theory is the AdS/CFT correspondence, to 
be discussed later. 

(iv) Assuming (local) thermal equilibrium, and hence the formation of a quark-gluon plasma 
(QGP), the question is whether this plasma is weakly or strongly coupled. The maximal temperature 
of this plasma, as estimated from the average energy density, should be around T ~ 500 600 MeV; 
so the respective coupling is moderately strong: a s ~ 0.3 -r- 0.4 or g ~ 1.5 -r- 2. The thermodynamic 
properties (like pressure or energy density) of a QGP at global thermal equilibrium within this range of 
temperatures are by now well known from numerical calculations on a lattice and can serve as a baseline 
of comparison for various effective theories. If the coupling is weak, one has to use the Hard Thermal 
Loop effective theory (HTL), a version of the kinetic theory which describes the long-range (or 'soft') 
excitations of the QGP. This effective theory lies at the basis of a physical picture of the QGP as a gas of 
weakly-coupled quasi-particles — quarks or gluons with temperature-dependent effective masses and 
couplings. Using this picture as a guideline for reorganizations of the perturbation theory, one has been 
able to reproduce the lattice data quite well. Thus, the thermodynamics appears to be consistent with 
a weak-coupling picture for the QGP, although this picture is considerably more complicated than that 
emerging from naive perturbation theory (the strict expansion in powers in g). Yet, this is not the end of 
the problem, as we shall see. 

(v) The QGP created in the intermediate stages of a HIC is certainly not in global thermal equi- 
librium, but only in a local one: it keeps expanding. Under very general assumptions, the effective 
theory describing this flow is hydrodynamics. The corresponding equations of motion are simply the 
conservation laws for energy, momentum, and other conserved quantities (like the electric charge or the 
baryonic number), and as such they are universally valid. But these equations also involve 'parameters', 
like the viscosities, which describe dissipative phenomena occurring during the flow and which depend 
upon the specific microscopic dynamics. The values of these parameters are very different at weak vs. 
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strong coupling. A meaningful way to characterize the strength of dissipation is via the dimensionless 
viscosity-over-entropy-density ratio 77/5. (This ratio is dimensionless when using natural units; other- 
wise it has the dimension of h.) Remarkably, the elliptic flow data at RHIC and the LHC to be later 
discussed suggest a very small value for this ratio, which is inconsistent with the present calculations 
at weak coupling, based on kinetic theory. On the other hand, such a small ratio is naturally emerg- 
ing at strong coupling, as shown by calculations within the AdS/CFT correspondence. The smallness of 
r]/s represents so far the strongest argument in favour of a strongly coupled quark-gluon plasma (sQGP). 
This may look contradictory with the previous conclusions drawn from thermodynamics. But one should 
remember that the QCD coupling depends upon the relevant space-time scale and that hydrodynamics 
refers to the long-range behaviour of the fluid, as encoded in its softest modes. By contrast, thermody- 
namics is rather controlled by the hardest modes — those with typical energies and momenta of the order 
of the (local) temperature. So, it is not inconceivable that a same system look effectively weakly coupled 
for some phenomena and strongly coupled for some others. 

(vi) Another strategy for studying the hadronic matter produced in a HIC refers to the use of hard 
probes. These are particles with large transverse energies (say, E± > 20 GeV at the LHC), which are 
produced in the very first instants of a collision and then cross the QCD matter liberated at later stages 
along their way towards the detector. Some of these particles, like the (direct) photons and the dilepton 
pairs, do not interact with this matter and hence can be used a baseline for comparaisons. But other 
particles, like quarks, gluons, and the jets initiated by them, do interact, and by measuring the effects of 
these interactions — say, in terms of energy loss, or the suppression of multi-particle correlations — one 
can infer informations about the properties of the matter they crossed. The RHIC data have demonstrated 
that semi-hard partons can lose a substantial fraction of the transverse energy via interactions in the 
medium ('jet quenching'), thus suggesting that these interactions can be quite strong. These results have 
been confirmed at the LHC, which moreover found that even very hard jets (E± > 100 GeV) can be 
strongly influenced by the medium, in the sense that they get strongly defocused : the energy distribution 
in the polar angle with respect to the jet axis becomes much wider after having crossed the medium. 
This is visible for the photon-hadron di-jet event in the right panel of Fig. [2] : the photon and the parton 
which has initiated the hadronic jet have been created by a hard scattering, so they must have been 
balanced in transverse momentum at the time of their creation. Yet, the central peak in the hadronic 
jet, which represents the final 'jet' according to the conventional definition, carries much less energy 
than the photon jet. This is interpreted as the result of energy transfer to large polar angles (outside 
the conventional 'jet' definition) via in-medium interactions. In order to study such interactions, in 
particular high-density effects like multiple scattering and coherence, it is again useful to build effective 
theories. In that case too, it is not so clear whether the physics is controlled by mostly weak coupling, 
or a mostly strong one. (By itself, the jet is hard, but its coupling to the relatively soft constituents of 
the medium may be still governed by a moderately strong coupling.) In fact, the unexpectedly strong jet 
quenching observed at RHIC is sometimes interpreted as another evidence for strong coupling behaviour. 
Moderately strong coupling turns out to be the most difficult situation to deal with, so in these lectures we 
shall rather describe the effective theories proposed in the limiting situations of weak and, respectively, 
strong coupling. Within perturbative QCD, this is known as medium-induced gluon radiation. At strong 
coupling, it again relies on AdS/CFT. 

3 The Color Glass Condensate 

This chapter is devoted to the early stages of an ultrarelativistic heavy ion collision (HIC), that is, the 
wavefunctions of the energetic nuclei prior to the collision and the partonic matter liberated by the col- 
lision. As already mentioned, these early stages are the realm of high-density, coherent, forms of QCD 
matter, characterized by high gluon occupation numbers. Such forms of matter can be described in terms 
of strong, semi-classical, colour fields. In what follows, we shall explain this theoretical description, 
starting with the perhaps more familiar parton picture of QCD scattering at high energy. 
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Fig. 3: Left: a cartoon of the proton structure in its rest frame. Right: cartoon of one saturation disk in the infinite 
momentum frame (this will be discussed in Section \377\ ). 

3.1 The QCD parton picture 

The microscopic structure of a hadron depends upon the resolution scales which are used to probe it, 
that is, upon the kinematics of the scattering process. It furthermore depends upon the Lorentz frame in 
which the hadron is seen: unlike physical observables, like cross-sections, which are boost invariant, the 
physical interpretation of these observables in terms of partons depends upon the choice of a frame. This 
is best appreciated by first looking at a hadron (say, a proton) in its rest frame (RF), where the proton 
4-momentum reads P^ = (M, 0, 0, 0). The proton has the quantum numbers of a system of three quarks 
— the 'valence quarks' — which are bound by confinement in a colour singlet state. But this binding 
proceeds via the exchange of gluons, which in turn can generate additional quark-antiquark pairs (see 
Fig. [3]). All these partons are 'virtual', meaning that they keep appearing and disappearing, and have 
typical energies and momenta of order Aqcd, since this is the scale where the QCD coupling becomes 
of 0(1) and thus the binding is most efficient. Clearly, such fluctuations are non-perturbative. Aqcd is 
also the typical scale for vacuum fluctuations, like a quark-antiquark pair pumping up from the vacuum 
and then being reabsorbed. By the uncertainty principle, such fluctuations have lifetimes and sizes of 
order 1/Aqcd, of the same order as the proton size itself. Under these conditions, it makes no sense to 
speak about 'hadronic substructure' : the hadronic fluctuations are ephemeral, delocalized over the whole 
proton volume, and cannot be distinguished from the vacuum fluctuations having the same kinematics 
and quantum numbers. 



Fig. 4: A hadronic fluctuation in the hadron rest frame (left) and in the infinite momentum frame (right). 

However, the situation changes if one observes the same hadron in a frame which is boosted by a 
large Lorentz factor 7 ^> 1 w.r.t. the rest frame. Then the hadron 4-momentum reads P M = (£", 0, 0, P) 
with E = \/P 2 + M 2 ~ P. (We have chosen the boost along the z axis and denoted P z = P.) In 
this boosted frame, conventionally referred to as the infinite momentum frame (IMF), the lifetime of the 
hadronic fluctuations is enhanced by Lorentz time dilation (see Fig. [4]), 





P = yP »M 
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Fig. 5: Left: the DIS process. Right: the absorption of the virtual photon by a quark. 

so these fluctuations are now well separated from the those of the vacuum (which have a lifetime 
~ 1/Aqcd in any frame, since the vacuum is boost invariant). The lifetime ^ is much larger than 
the duration of a typical collision process (see below); so, for the purpose of scattering, the hadronic 
fluctuations can be viewed d&free, independent quanta. These quanta are the partons (a term coined by 
Feynman). It then becomes possible to factorize the cross-section (say, for a hadron-hadron collision) 
into the product of parton distribution functions (one for each hadron partaking in the collision), which 
describe the probability to find a parton with a given kinematics inside the hadronic wavefunction, and 
partonic cross-sections, which, as their name indicates, describe the collision between subsets of par- 
tons from the target and the projectile, respectively. If the momentum transferred in the collision is hard 
enough, the partonic cross-sections are computable in perturbation theory. The parton distributions are 
a priori non-perturbative, as they encode the information about the binding of the partons within the 
hadron. Yet, there is much that can be said about them within perturbation theory, as we shall explain. 
To that aim, one needs to better appreciate the role played by the resolution of a scattering process. In 
turn, this can be best explained on the example of a simpler process: the electron-proton deep inelastic 
scattering (DIS). 

The DIS process is illustrated in Fig. [5] (left): an electron with 4-momentum scatters off the 
proton by exchanging a virtual photon (7*) with 4-momentum and emerges after scattering with 
4-momentum = — q^. The exchanged photon is space-like : 

q 2 = (£- I'f = -21 • g = -2E £ E £ f(l - cos 0a') = ~Q 2 with Q 2 > 0, (3) 

with Ei = \£\, E11 = \£ f \, and Qn> = /.{£,£'). The positive quantity Q 2 is referred to as the 'virtuality'. 
The deeply inelastic regime corresponds to Q 2 » M 2 , since in that case the proton is generally broken 
by the scattering and its remnants emerge as a collection of other hadrons (denoted by X in Fig. [5]). The 
(inclusive) DIS cross-section involves the sum over all the possible proton final states X for a given £' '. 

A space-like probe is very useful since it is well localized in space and time and thus provides 
a snapshot of the hadron substructure on controlled, transverse and longitudinal, scales, as fixed by the 
kinematics. Specifically, we shall argue that, when the scattering is analyzed in the proton IMF, the 
virtual photon measures partons which are localized in the transverse plane within an area S ~ 1/Q 2 
and which carry a longitudinal momentum k z = xP, where x is the Bjorken variable : 

Q 2 Q 2 
X = 2(P • q) ~ s + Q 2 - M 2 ' (4) 

where s = (P + q) 2 is the invariant energy squared of the photon+proton system. That is, the two 
kinematical invariants Q 2 and x, which are fixed by the kinematics of the initial state (£, P) and of the 
scattered electron (l f ), completely determine the transverse size (~ l/Q) and the longitudinal momentum 
fraction (x) of the parton that was involved in the scattering. This parton is necessarily a quark (or 
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antiquark), since the photon does not couple directly to gluons. But the DIS cross-section allows us to 
indirectly deduce also the gluon distribution, as we shall see. 

As a first step in our argument, consider a quark excitation of the hadron, viewed in the IMF. 
This quark is a virtual fluctuation which has been boosted together with the proton, so its virtuality and 
its transverse momentum are both small as compared to its longitudinal momentum k z = £P. (We 
temporarily denote with £ the fraction of the proton longitudinal momentum which is carried by the 
quark.) So, for most purposes, one can treat the quark as a nearly on-shell excitation with 4-momentum 
as k^ ~ = (£P,0 ? 0,£P) and k 2 = k^ ~ 0. (More precisely, k 2 ~ k\ - A^ CD .) Such an 
excitation has a relatively large lifetime, which can be estimated as in Eq. ([2]) : 

2k z 2£P 

At fluct 7A£rf - T2" = "T2- j (5) 

where A£rf ^ is the lifetime of the fluctuation in the hadron rest frame and 7 = k z /k± is the 

boost factor from the RF to the IMF. 

Consider now the absorption of the virtual photon by the quark, cf . Fig. [5] right. The quark is 
liberated by this collision, meaning that it is put on shell; so we can write 

(k + qf = => -Q 2 + 2£P.q = £ = = x, (6) 

where we have also used k 2 « 0, as discussed before. We see that the collision identifies the longitudinal 
momentum fraction £ of the participating quark with the Bjorken-x kinematical variable, as anticipated. 
From now on, we shall use the notation x for both quantities. 

To also clarify the transverse resolution of the virtual photon, we first need an estimate for the 
collision time. This is the typical duration of the partonic process q + 7* — >> q (cf. Fig. fright) and 
is given by the uncertainty principle: At co n ~ 1/ AE, where AE = qo + |fc + q| — |fc| is the energy 
difference at the photon emission vertex. To estimate AE, it is convenient to choose a space-like photon 
with zero energy and only transverse momentum: q^ = (0, q±, 0). Then 



/ £7 2xP 

AE=\k + q\- \k\ = ^(xP) 2 + q\ - xP - ^ =^ At co n ~ (7) 

(Note that Q 2 = q\ for the virtual photon at hand.) In order to be 'found' by the photon, a quark 
excitation must have a lifetime larger than this collision time: 

2xP 2xP 
At fluct - -p- > At co n ~ -^2- fe± < Q 2 • (8) 

Hence, the virtual photon can discriminate only those partons having transverse momenta smaller than 
its virtuality Q. By the uncertainty principle, such partons are localized within a transverse area ~ l/Q 2 , 
as anticipated after Eq. Q. 

The previous considerations motivate the following formula for the DIS cross-section : 

^p(x,Q 2 ) = ^^F 2 (x,Q 2 ), (9) 

where the first factor in the r.h.s. is the elementary cross-section for the photon absorbtion by a quark (or 
an antiquark), whereas the second factor — the structure function F2(x,Q 2 ) — is the sum of the quark 
and antiquark distribution functions, weighted by the respective electric charges squared 

F 2 (x,Q 2 ) = J2 e } [xqf(x,Q 2 ) + xq f (x,Q 2 )] , 
/ 
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That is, qf(x, Q 2 )dx is the number of quarks of flavor / with longitudinal momentum fraction between 
x and x + dx and which occupy a transverse area l/Q 2 . 

One may naively think that the condition k\ < Q 2 is trivially satisfied, since the partons confined 
inside the hadron have transverse momenta k± ~ Aqcd> whereas Q 2 » ^qcd by definition of 
DIS. If that were the case, the structure function ^(x, Q 2 ) would be independent of Q 2 — a property 
known as Bjorken scaling. However, the DIS data show that Bjorken scaling holds only approximately 
and only in a limited range of values for x, namely for x > 0.1. This can be understood as follows : the 
typical transverse momenta are ~ Aqcd only for the valence quarks and, more generally, for the partons 
with relatively large longitudinal momentum fractions. But virtual quanta with much larger values for 
k± can be generated via radiative processes like bremsstrahlung. Such quanta have very short lifetimes, 
but so long as k\ < Q 2 , they can still contribute to DIS. Also, they generally have small values of x, as 
they share all together the longitudinal momenta of their parents partons. Hence, we expect the parton 
evolution via bremsstrahlung to lead to an increase in the parton distributions at large values of k± and 
small values of x. This evolution is responsible for the violations of the Bjorken scaling seen in the data 
and, more generally, for the DGLAP evolution |I][2j of the parton distribution functions with increasing 
Q 2 . It is furthermore responsible for the rapid growth in the gluon distribution with decreasing x and the 



formation of a colour glass condensate at high energy. This will be further discussed in Section 3.3 



3.2 Particle production at the LHC: why small x ? 

Before we turn to a discussion of parton evolution, let us explain here why we shall be mostly interested 
in partons with small longitudinal momentum fractions x <C 1. As it should be clear from Eq. ([4]), 
small values of x correspond to the high-energy regime at s » Q 2 . The conceptual importance of this 
regime will be explained later, but for the time being let us discuss it from the experimental point of view. 
Very small values of x, as low as x = 10 -5 10 -4 , have been already reached in the e+p collisions 
at HERA, but in that context they were associated (because of the experimental constraints) with rather 
small values of the transferred momentum Q 2 . Namely, the HERA data at x < 10 -4 correspond to 
values Q 2 < 1 GeV 2 which are only marginally under control in perturbation theory. Because of that, 
the DIS data at HERA remained inconclusive for a check of our theoretical understanding of the physics 
at small x. 

But the situation has changed with the advent of the new hadron-hadron colliders, RHIC and, 
especially, the LHC. Given the much higher available energies, the bulk of the particle production (with 
semi-hard transverse momenta) in these experiments is controlled by partons with x < 10 -3 . Moreover, 
for special kinematical conditions to be shortly specified, one can probe values as low as x ~ 10 -6 with 
truly hard momentum transferts, such as Q 2 = 10 GeV 2 . 

To describe the kinematics of particle production, it is useful to introduce a new kinematical vari- 
able, the rapidity y, which is an alternative for the longitudinal momentum. For an on-shell particle with 
4-momentum p^ = (E,p±,p z ), the rapidity is defined as 

y = - In — — ^ z =^ E — m± cosh y. p z — m± sinh y , (11) 
2 E — p z 

where m± = ^Jm 2 + p\ is the 'transverse mass' and E 2 = m 2 ^ + p 2 . Note that y is positive for a 
'right-mover' (p z > 0) and negative for a 'left-mover' (p z < 0). In fact, one has v z = p z j E = tanhy, 
so y is simply related to the longitudinal boost factor: 7 = coshy. A similar quantity which is perhaps 
more useful in the experiments (since easier to measure) is the pseudo-rapidity 

rj = - In - In tan - , (12) 

2 p-Pz 2 
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Fig. 6: Left: two particle production in hadron-hadron collision. Right: the kinematical domain accessible at the 
LHC, as compared to RHIC and the SPS; the small-x region is highlighted. 



where p = \p\is the magnitude of the 3-momentum vector. As shown by the second equality above, 77 is 
directly related to the polar angle 9 made by the particle with the longitudinal axis (cos 9 = p z /p)> For 
massless particles or for ultrarelativistic ones (whose masses can generally be ignored), the two rapidities 
coincide with each other, as manifest by comparing Eqs. ( [TT] ) and ( [T2| ). 

Consider now the process illustrated in Fig. [6] (left), i.e. the production of a pair of particles in 
a partonic subcollision of a hadron-hadron scattering. In the center-of-mass (COM) frame, the two 
partons partaking in the collision have 4-momenta = XiPj 1 + where i = 1, 2, P± = (P, 0, 0, P), 
p% = (p, 0, 0, -P), and k? ± = (0, k i± , 0). Notice that P = y/s/2. The two outgoing particles will 
be characterized by the respective transverse momenta, p a ± and pb±, and rapidities, y a and y&. Energy- 
momentum conservation implies p a ± + p^± = k\± + &2± and 

Xl = ^ e y + ^ e"> , x 2 = ^ e~y* + ^ e"^ . (13) 

For particle production at RHIC or the LHC, the average transverse momentum of a hadron in the final 
state is below 1 GeV ; moreover, 99% of the 'multiplicity' (i.e. of the total number of produced hadrons) 
has p± < 2 GeV (see Fig. [7]). For p a ^_\_ = 1 GeV and central rapidities y a ^ ~ 0, Eq. ( [T3| ) implies 

Xi ~ 1CT 2 at RHIC (y/s = 200 GeV), x { ~ 4 x 10" 4 at the LHC (y/s = 2.76 TeV), (14) 

where in the case of the LHC we have chosen the maximal COM energy per nucleon pair that has been 
reached so far in Pb+Pb collisions. Thus, the bulk of the particle production is initiated by partons 
carrying small values of x, as anticipated. Moreover, one of these values (x\ or xi) can be made much 
smaller by studying particle production at either forward, or backward, rapidities. The rapidities are 
'forward' when both y a and are positive and relatively large, that is, the final particles propagate 
essentially along the same direction as the original hadron T ; then their production probes very small 
values of X2 in the wavefunction of the hadron '2' and comparatively large values of x\. At the LHC, 
one can probe values as small as X2 ~ 10 -6 for p± ~ 10 GeV, as indicated in the r.h.s. of Fig. [6] 

We finally discuss the cross-section for the production of a pair of hadrons. When the transverse 
momenta p a ^ ± are large enough, one can ignore the 'intrinsic' transverse momenta j_ °f the colliding 
partons. Then the transverse momentum conservation p a ± + p^± ~ implies that the outgoing particles 
propagate back-to-back in the transverse plane, i.e. they make an azimuthal angle Acf) ~ tt. The 
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Fig. 7: Transverse momentum dependence of the single inclusive particle production in hadron-hadron collisions. 
Left: Pb+Pb collisions at the LHC. Right: p+p collisions at the LHC and p+p collisions at SPS and CDF. 



associated cross-section admits the following collinear factorization, analogous to Eq. fL0| ) for DIS 

= ^ xifi{xi,ii 2 ) x 2 fj(x 2l /i 2 ) , (15) 



dp 2 L dyidy 2 ^ * JJX ' dp 



where xfi(x, /r ) are parton distributions for all species of partons (i = q, q, g), fi 2 is the factorization 
scale, and da^ /dp^ is the cross-section for the (relatively hard) partonic process i+j — » a+6. Leading- 
order perturbative QCD yields da/dp^ oc o? s jp\_ at high energy. So, if one tries to compute the total 
multiplicity by integrating over all values of p\ (say, for y a ^ Vb ^ 0), then one faces a quadratic 
infrared divergence from the limit p\ — )► 0. One may think that this divergence is cut off at p± ~ Aqcd, 
since this is the typical value expected for the intrinsic momenta ki^±> But then one would conclude 
that the bulk of the particle production, even at very high energy, is concentrated at very soft transverse 
momenta, of the order of the confinement scale Aqcd- Moreover, the average p± would be independent 
of the energy (since of O(Aqcd))- These conclusions are however contradicted by the data in Fig. [7j 
which rather show that (p±) ^0.5 GeV is about 2 to 3 times larger than Aqcd at the LHC energies 
and, remarkably, it clearly rises with the COM energy E = y/s. This conflict between the data and 
the prediction ( [T5] ) of collinear factorization clearly shows that the latter cannot be extrapolated down 
to lower values for p±, say of order 1 GeV. The proper way to describe this semi-hard region within 
(perturbative) QCD will be explained in the next subsection. The main outcome of that analysis will be 
to introduce a new infrared cutoff in the problem, which is dynamically generated — via gluon evolution 
with decreasing x — and rises as a power of the energy. This is the saturation momentum. 

3.3 Gluon evolution at small x 

In perturbative QCD, parton evolution proceeds via bremsstrahlung, which favors the emission of soft 
and collinear gluons, i.e. gluons which carry only a small longitudinal momentum fraction x«l and 
a relatively small transverse momentum k±. Fig. [8] illustrates one elementary step in this evolution: 
the emission of a gluon which carries a fraction x = k z /p z of the longitudinal momentum of its parent 
parton (quark or gluon). For x«l and to lowest order in a s , the differential probability for this emission 
(obtained as the modulus squared of the amplitude represented in Fig. [8]) reads 

dP B rem - Cr ^~ , (16) 

7T Z K\ X 
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J&flflQflflQflQflfl. / k z = xp z 5 k ± 

Fig. 8: Gluon bremsstrahlung out of a parent quark to lowest order in pQCD. 

where Cr is the SU(iV c ) Casimir in the colour representation of the emitter: Ca = ^V c for a gluon and 
C F = (N 2 — 1)/2N C for a quark. (N c is the number of colours, which is equal to 3 in real QCD, but it 
is often kept as a free parameter in theoretical studies, because many calculations simplify in the formal 
limit N c » 1. The results obtained in this limit provide insightful, qualitative and semi-quantitative 
informations about real QCD.) Eq. ( fT6] ) exhibits the collinear (k± —> 0) and soft (x —> 0) singularities 
mentioned above, which result in the enhancement of gluon emission at small k± and/or x. If the emitted 
parton with small x were a quark instead a gluon, there would be no small x enhancement, only the 
collinear one. This asymmetry, due to the spin-1 nature of the gluon, has the remarkable consequence 
that the small-x part of the wavefunction of any hadron is built mostly with gluons. 



As manifest on Eq. (16), parton branching is suppressed by a power of a s {k 2 A ), which is small 
when k± » Aqcd- But this suppression can be compensated by the large phase-space available for the 
emission, which equals ln(Q 2 /AQ CD ) for the emission of a parton (quark or gluon) with transverse mo- 
mentum k± <C Q and, respectively, ln(l/x) for that of a gluon with longitudinal momentum fraction £ 
within the range i<^<1. Hence, for large Q 2 » Aq CD and/or small xCl, such radiative processes 
are not suppressed anymore and must be resummed to all orders. Depending upon the relevant values 
of Q 2 and x, one can write down evolution equations which resum either powers of a s \n(Q 2 / Aq CD ), 
or of a s \n(l/x), to all orders. The coefficients in these equations represent the elementary splitting 
probability and can be computed as power series in a s , starting with the leading-order result in Eq. fL6| ). 

The evolution with increasing Q 2 is described by the DGLAP equation (from Dokshitzer, Gribov, 
Lipatov, Altarelli and Parisi) |TJ|2j. This evolution mixes quarks and gluons (see Fig. [9ja), which in 
particular allows us to reconstruct the gluon distribution from the experimental results for F^. The small- 
x evolution, on the other hand, involves only gluons and corresponds to resumming ladder diagrams like 
those in Fig.[9jb in which successive gluons are strongly ordered in x (see below). Both evolutions lead 
to an increase in the number of partons at small values of x (and a decrease at large values x > 0.1), but 
the physical consequences are very different in the two cases: 

(i) When increasing Q 2 , one emits partons which occupy a smaller transverse area ~ l/Q 2 , as 



shown in Fig. 11 (right). The decrease in the area of the individual partons is much stronger than the 
corresponding increase in their number. Accordingly, the occupation number in the transverse plane 
decreases with increasing Q 2 , meaning that the partonic system becomes more and more dilute. Accord- 
ingly, the partons may be viewed as independent. This observation lies at the basis of the conventional 
parton picture, which applies for sufficiently high Q 2 (at a given value of x). 

The parton occupation number mentioned above yields the proper measure of the parton density in 
the hadron. It can be estimated as [the number of partons with a given value ofx] x [the area occupied 
by one parton] divided by [the transverse area of the hadron], that is (for gluons, for definiteness), 



"<*• = i?^ ■ w,th s{x - Q ] = -I? iQ ) ~ J d k± drffc • <17) 

where R is the hadron radius in its rest frame (so its transverse area is ~ ttR 2 in any frame). The 
numerator in the above definition of the occupation number, that is 

/ dN q , dN q AN q 

xg (x,Q 2 = x— ± = k z — ~ — — 9 —, (18) 
dx ak z Az Ak z 
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Fig. 9: Parton evolution in perturbative QCD. Theparton cascade on the right involves only gluons (at intermediate 
stages) and is a part of the BFKL resummation at small x (see the discussion around Eq. (22fy). 



is known as the gluon distribution. The last estimate above follows from the uncertainty principle: 
partons with longitudinal momentum k z = xP are delocalized in z over a distance Az ~ l/k z . Hence, 
the gluon distribution yields the number of gluons per unit of longitudinal phase-space, which is indeed 
the right quantity for computing the occupation number. Note that gluons with x <C 1 extends in z 
over a distance Az ~ 1/xP which is much larger than the Lorentz contracted width of the hadron, 
R/j ~ 1/P. This shows that the image of an energetic hadron as a 'pancake', that would be strictly 
correct if the hadron was a classical object, is in reality a bit naive: it applies for the valence quarks with 
x ~ 0(1) (which carry most of the total energy), but not also for the small-x partons (which are the 
most numerous, as we shall shortly see). 

(ii) When decreasing x at a fixed Q 2 , one emits mostly gluons which have smaller longitudinal 
momentum fractions, but which occupy, roughly, the same transverse area as their parent gluons (see 



Fig. [TT] right). Then the gluon occupation number, Eq. ( |T7| ), increases, showing that the gluonic system 
evolves towards increasing density. As we shall see, this evolution is quite fast and eventually leads to a 
breakdown of the picture of independent partons. 

In order to describe the small-x evolution, let us start with the gluon distribution generated by 
a single valence quark. This can be inferred from the bremsstrahlung law in Eq. fl6| ) (the emission 
probability is the same as the number of emitted gluons) and reads 

dN g 2 _ a s C F [ Q2 dk\ a s C F n ( Q 2 



dx n A2 CD k\ it V A QCd/ 

where we have ignored the running of the coupling — formally, we are working to leading order (LO) in 
pQCD where the coupling can be treated as fixed — and the 'infrared' cutoff Aqcd has been introduced 
as a crude way to account for confinement: when confined inside a hadron, a parton has a minimum 
virtuality of (9(Aq CD ). In Eq. (19) it is understood that x <C 1. In turn, the soft gluon emitted by 
the valence quark can radiate an even softer gluon, which can radiate again and again, as illustrated in 



figure [TOj Each emission is formally suppressed by a power of a s , but when the final value of x is tiny, 
the smallness of the coupling constant can be compensated by the large available phase-space, of order 
ln(l/x) per gluon emission. This evolution leads to an increase in the number of gluons with 

For a quantitative estimate, consider the first such correction, that is, the two-gluon diagram in 



Fig. 10 left: the region in phase-space where the longitudinal momentum fraction x\ of the intermediate 
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Fig. 10: Gluon cascades produced by the high-energy (BFKL) evolution of the proton wavefunction. 



gluon obeys x <C x\ <C 1 provides a contribution of relative order 

a s N c f 1 dxi 1 _ ot s N c 

/ = a s In - , «.s = • (20) 

7T J x X\ X TT 

When a s \n(l/x) ~ 1, this becomes of 0(1), meaning that this two-gluon diagram contributes on the 
same footing as the single gluon emission in Fig.[8j A similar conclusion holds for a diagram involving n 



intermediate gluons strongly ordered in x, cf. Fig. 10 right, which yields a relative contribution of order 



a ;f^f^...f> I^IV, <21) 

Jx x n Jx n x n—l J X2 X l Ti. \ X J 



When a s \n{l/x) > 1, the correct result for the gluon distribution at leading order is obtained by sum- 
ming contributions from all such ladders. As clear from Eq. ( [2T] ), this sum exponentiates, modifying the 



integrand of Eq. ( 19 ) into 



W^° Y > Y ^~> (22) 
dxdk^ TT fcj_ x 

where uj is a number of order unity which cannot be determined via such simple arguments. The variable 
Y is the rapidity difference between the final gluon and the original valence quark and it is often simply 



referred to as 'the rapidity'. The quantity in the l.h.s. of Eq. ( |22| ) is the number of gluons per unit rapidity 
and with a given value k± for the transverse momentum, a.k.a. the unintegrated gluon distribution^ 

To go beyond this simple power counting argument, one must treat more accurately the kinematics 
of the ladder diagrams and include the associated virtual corrections. The result is the BFKL equation 
(from Balitsky, Fadin, Kuraev, and Lipatov) [ 3 ] for the evolution of the unintegrated gluon distribution 
with Y. The solution of this equation, which resums perturbative corrections (a s Y) n to all orders, 



confirms the exponential increase in Eq. ( |22| ), albeit with a fc^-dependent exponent and modifications to 
the k J 2 -spectrum of the emitted gluons. 

An important property of the BFKL ladder is its coherence in time : the lifetime of a parton being 
proportional to its value of x, At ~ 2k z /k\ oc x, cf. Eq. ([5]), the 'slow' gluons at the lower end of the 
cascade have a much shorter lifetime than the preceding 'fast' gluons. Therefore, for the purposes of 
small-x dynamics, fast gluons with x' ^> x act as frozen colour sources emitting gluons at the scale x. 
Because these sources may overlap in the transverse plane, their colour charges add coherently, giving 
rise to a large colour charge density. The average colour charge density is zero by gauge symmetry 
but fluctuations in the colour charge density — as measured in particular by the unintegrated gluon 



distribution — are nonzero and increase rapidly with 1/x, cf. Eq. (22). 



2 The occupation number {XT} is more correctly defined as the unintegrated gluon distribution per unit transverse area: 
n(Y, k±) = dN 9 / (dYd 2 k±d b±) where b± (the 'impact parameter') is the transverse position of a gluon with respect to the 
center of the hadron. 
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Fig. 11: Left: the 1 / x-evolution of the gluon, sea quark, and valence quark distributions for Q 2 = 10 GeV 2 , as 
measured at HERA (combined HI and ZEUS analysis Note that the gluon and sea quark distributions have 
been reduced by a factor of 20 to fit inside the figure. Right: the 'phase-diagram' for parton evolution in QCD; 
each coloured blob represents a parton with transverse area Ax±_ ~ 1/Q 2 and longitudinal momentum k z = xP. 
The straight line In Q 2 (x) = XY is the saturation line, cf. Eq. (25), which separates the dense and dilute regimes. 



This growth is indeed seen in the data: e.g., the HERA data for DIS confirm that the proton 



wavefunction at x < 0.01 is totally dominated by gluons (see Fig. [TT] left). However, on physical 
grounds, such a rapid increase in the gluon distribution cannot go on for ever (that is, down to arbitrarily 
small values of x). Indeed, the BFKL equation is linear — it assumes that the radiated gluons do not 
interact with each other, like in the conventional parton picture. While such an assumption is perfectly 
legitimate in the context of the Q 2 -evolution, which proceeds towards increasing diluteness, it eventually 
breaks down in the context of the y-evolution, which leads to a larger and larger gluon density. As long 



as the gluon occupation number ( [17] ) is small, n <C 1, the system is dilute and the mutual interactions 
of the gluons are negligible. When n ~ 0(1), the gluons start overlapping, but their interactions are 
still weak, since suppressed by a s <C 1. The effect of these interactions becomes of order one only 
when n is as large as n ~ 0(l/a s ). When this happens, non-linear effects (to be shortly described) 
become important and stop the further growth of the gluon distribution. This phenomenon is known as 
gluon saturation pHT). An important consequence of it is to introduce a new transverse-momentum 



scale in the problem, the saturation momentum Q s (x), which is determined by Eq. ( 17 ) together with the 
condition that n ~ l/a s : 

n(x, Q 2 = Ql{x)) ~ i- =► Ql(x) c a /^f . (23) 

Except for the factor a s , the r.h.s. of Eq. ( [23] ) is recognized as the density of gluons per unit transverse 
area, for gluons localized within an area S ~ 1/Q 2 s {x) set by the saturation scale. Gluons with k± < 
Q s (x) are at saturation: the corresponding occupation numbers are large, n ~ l/a s , but do not grow 
anymore when further decreasing x. Gluons with k± » Q s {x) are still in a dilute regime: the occupation 
numbers are relatively small n l/a s , but rapidly increasing with 1/x via the BFKL evolution. The 
separation between the saturation (or dense, or CGC) regime and the dilute regime is provided by the 



saturation line in Fig. 1 1 right, to be further discussed below. 



The microscopic interpretation of Eq. ( |23| ) can be understood with reference to Fig. [12] (left) : 
gluons which have similar values of x (and hence overlap in the longitudinal direction) and which occupy 
a same area ~ l/Q 2 in the transverse plane can recombine with each other, with a cross-section o- gg ^ g ~ 
a s /Q 2 . After taking also this effect into account, the change in the gluon distribution in one step of the 
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Fig. 12: Left: gg —> g recombination process leading to saturation. Right: the saturation momentum Q 3 (x, A) as 
a function of the longitudinal momentum fraction x and of the atomic number A. 



small-x evolution (i.e. under a rapidity increment Y —> Y + dY) can be schematically written as 



d_ 
W 



xg(x,Q 2 ) = cja s xg(x,Q 2 



a* 



Q 2 R 2 



[xg(x, Q 



2M2 



(24) 



The overall factor of a s in the r.h.s. comes from the differential probability oc a s dY to emit one addi- 
tional gluon in this evolution step, cf. Eq. ( 16). The first term, linear in xg{x, Q 2 ), represents the BFKL 
evolution; by itself, this would lead to the exponential growth with Y shown in Eq. ( [22] ). The second 
term, quadratic in xg{x 1 Q 2 ), is the rate for recombination. This is formally suppressed by one factor a s , 
but it becomes as important as the first term when Q 2 is of the order of the saturation momentum Q 2 (x) 
introduced in Eq. ( [23] ). When that happens, the r.h.s. of Eq. ( [24] ) vanishes, and then the gluon distribu- 
tion stops growing with Y. The above argument, due to Gribov, Levin and Ryskin back in 1983 (5j, is 
a bit oversimplified (and the actual evolution equation is considerably more complicated than Eq. ( [24] ); 
see the review papers |8]-[T5| and the discussion in Section [374] below), but it has the merit to illustrate 
in a simple way the physical mechanism at work: the gluon occupation numbers saturate because the 
non-linear effects associated with the high gluon density compensate the bremsstrahlung processes. 



Remarkably, Eq. (23 ) implies that the saturation momentum increases with 1/x, since so does the 



gluon distribution for k± > Q s (x), cf. Eq. (22). So, for sufficiently small values of x (say, x < 10 in 



the case of a proton), one expects Q 2 s {x) » Aqcd- In that case, the (semi)hard scale Q s (x) supplants 
Aqcd as an infrared cutoff for the calculation of physical observables like the multiplicity (cf. the 
discussion at the end of Section [T2] ). This has the remarkable consequence that, for sufficiently high 
energy, the bulk of the particle production can be computed in perturbation theory. But the proper 
framework to perform this calculation is not standard pQCD as based on the collinear factorization, 
but the CGC effective theory which includes the non-linear physics of gluon saturation. This will be 
discussed in the next subsection. 

Gluon occupancy is further amplified if instead of a proton we consider a large nucleus with 
atomic number A » 1. The corresponding gluon distribution xgA(x, Q 2 ) scales like A, since gluons 
can be radiated by any of the 3 A valence quarks of the A nucleons. Since the nuclear radius scales like 
Ra ~ A 1 / 3 , Eq. implies that the gluon occupation number scales as A 1 / 3 . This factor is about 6 
for the Au and Pb nuclei respectively used at RHIC and the LHC. Thus, for a large nucleus, saturation 
effects become important at larger values of x than for a proton. This explains why ultrarelativistic heavy 
ion collisions represent a privileged playground for observing and studying the effects of saturation. 



Fig. 12 (right) summarizes our current expectations for the value and the variation of the saturation 
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Fig. 13: Charged particles multiplicity at central (pseudo)rapidity r] = as a function of the COM energy forp+p 
and nucleus-nucleus collisions. The data are consistent with a power-law increase with y/s (with an exponent 
which is slightly larger for A+A than forp+p), but appear to exclude a logarithmic law oc In s. 



momentum. The dependence upon x is by now known to next-to-leading-order (NLO) accuracy fT6 | 
— that is, by resumming radiative corrections a s [a s \n(l/x)] n to all orders together with non-linear 
effects. The result can be roughly expressed as 

Q 2 s (x,A) ~ Ql Al/3 (^Y > with A = 0.2 -0.3, (25) 

with the power A known as the saturation exponent. The overall scale Qq, which has the meaning of 
the proton saturation scale at the original value xq, is non-perturbative and cannot be computed within 
the CGC effective theory. (The latter governs only the evolution from xo down to x <C x$.) In practice, 
this is treated as a free parameter which is fitted from the data. The fits yield Qq ~ 0.5 GeV for 



xq = 10 2 . Fig.[12jshows that for x = 10 b (a typical value for forward particle production at the LHC), 



-5 



Q s — 1 GeV for the proton, while Q s ~ 3 GeV for the Pb nucleus. This difference is significant: while 
1 GeV is only marginally perturbative, 3 GeV is sufficiently 'hard' to allow for controlled perturbative 
calculations. This confirms the usefulness of HIC as a laboratory to study saturation. 

Before we conclude this subsection, let us notice some robust predictions of the saturation physics, 
which do not require a detailed theory and can be directly checked against the data. One of them refers 
to the energy-dependence of the average transverse momentum of the produced particles: as shown in 
Fig.[7j this grows like a power of E = y/s, with an exponent which is fitted from the data as 0.115. This 
is consistent with expectations based on gluon saturation fT7| . Indeed, prior to the collision, the gluon 



distribution inside the hadron wavefunction is peaked at k± ~ Q s (see Fig. [16] below and the related 
discussion) and these gluons are then released in the final state. We thus expect the average p± of the 
produced hadrons to scale like Q s (x) evaluated at the appropriate value of x, that is, x = Pa/E (cf. 
Eq. ( 13 )). This argument implies (p±) oc E x / 2 , which is indeed consistent with the data in Fig. [7] (right) 



together with the estimate in ( [25] ) for the saturation exponent. Another prediction of this kind refers to 
the particle multiplicity in the final state d7V/d?7, say, at central (pseudo)rapidity 77 = 0. By the above 
argument, this is dominated by gluons with Q 2 ~ Ql(x) and hence it is proportional to the respective 



gluon distribution, that is, to Q s (x) itself (cf. Eq. (23)) : d7V/d?7 Ql(E) ~ E • Once again, this 



appears to be consistent with the data for both p+p and A+A collisions, as shown in Fig. 13 



3.4 The CGC effective theory 

The partonic form of matter made with the saturated gluons is known as the colour glass condensate 
~T5L 
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- This is coloured since gluons carry the 'colour' charge of the non-Abelian group SU(3). 

- It is a glass because of the separation in time scales, due to Lorentz time dilation, between the 
'slow' gluons at small x and their 'fast' sources at larger x. The sources appear as 'frozen' over 
the characteristic time scales for the dynamics at small x, but they can vary over much larger time 
scales, as set by their own, comparatively large, longitudinal momenta. A system which behaves 
as a solid on short time scales and as a fluid on much longer ones, is a glass. 

- It is a condensate because the saturated gluons and their sources have high occupation numbers 



n(x, k±) ~ l/a s and their colour charges add coherently to each other, as explained in Section|33 
in relation with the BFKL ladder. A coherent quantum state with high occupancy can be in a 
first approximation described as a classical field (here, a colour field), which is the most generic 
example of a condensate. 

Because of its high density, the CGC is weakly coupled and thus it can be studied within perturbative 
QCD. This is strictly correct for sufficiently small values of x, such that Q 2 (x) » Aq CD and hence 
a s (Q 2 ) <C 1, but it remains marginally true for the phenomenology at RHIC and, especially, the LHC, 
where the saturation momentum is semi-hard, cf. Fig. [12] (right). Based on that, an effective theory has 
been explicitly constructed, which resums an infinite series of Feynman graphs of the ordinary pertur- 
bation theory — those which are enhanced by either the large logarithm ln(l/x), or by the high gluon 
density. This theory governs the dynamics of the gluons with a given, small, value of x, while the gluons 
at larger values x r » x have been 'integrated out' in perturbation theory. In order to describe its mathe- 
matical structure, it is useful to recall that the gluon field in QCD is represented by a non-Abelian vector 
potential A%(x) where the upper index \i refers to the 4 Minkowski coordinates and the subscript a is a 
colour index in the adjoint representation of SU(iV c ) and can take N 2 — 1 = 8 values. 

The CGC effective theory may be viewed as a non-linear generalization of the BFKL evolution, 
but in fact it is much more complex than just a non-linear evolution equation (say, like that in Eq. ([24])). 
The BFKL equation applies to the unintegrated gluon distribution (or occupation number), which is a 
Fourier transform of the 2-point functioij^] (A l a (x)A l a (y)) of the colour fields within the hadron (The 
average refers to the hadron wavefunction and the upper index i with i = 1, 2 indicates the transverse 
directions.) This quantity offers more information than the standard parton distributions like xg(x, Q 2 ) 
— it also describes the distribution of gluons in transverse momentum, and not only in x — , but it 
still does not probe many-body correlations in the gluon distribution, as the higher n-point functions 
with n > 4 would do. The restriction to the 2-point function is justified so long as the system is 
dilute and gluons do not interact with each other. But this cannot encode the non-linear physics of 
saturation, which is sensitive to higher n-point functions and hence to correlations. In fact, to correctly 
describe gluon saturation, one needs to control n-point functions with arbitrarily high n. This can be 
understood as follows: the fact that the occupation numbers are n ~ 0(l/a s ) at saturation, means 
that the colour field strengths are as large as A l a ~ G(l/g), and then there is no penalty for inserting 
arbitrary powers of A l a . Indeed, any such an insertion is accompanied by a factor of g. (Recall that 
interactions in QCD enter via the covariant derivative = <9 M — igA^.) So, the CGC effective theory 
is truly an infinite hierarchy of coupled evolution equations describing the simultaneous evolution of all 
the n-point functions. Remarkably enough, this hierarchy can be summarized into a single, functional, 
evolution equation for the CGC weight function — a functional generalization of the 'unintegrated' gluon 
distribution that will be shortly discussed. 

The key ingredient for such an economical description is the proper choice of the relevant degrees 
of freedom: as already mentioned, the small-x gluons with high occupation numbers n ~ 1 /a s can be 
treated semi-classically to leading order in a s — that is, they can be described as classical colour fields 
Aa{x) radiated by colour sources representing the faster gluons with x r » x. This distinction between 
'classical fields' (= the small-x gluons for which the effective theory is built) and their 'sources' (= 



3 See Eq. d35|) for a more precise definition of the unintegrated gluon distribution in the presence of non-linear effects. 
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Fig. 14: Schematic representation of the d.o.f. involved in the CGC effective theory and of the quantum evolution 
which is taken into account in this theory. A newly emitted gluon with a small longitudinal momentum fraction 
x <C 1 rescatter off the gluon held [p] created in the previous steps by the gluons with larger values x' ^> x, 
effectively represented by their global colour charge density p a . 



the large-x gluons which are integrated out in the construction of the effective theory) is illustrated in 



Fig. 14 The effective theory based on this separation is valid to LO in a s , but to all orders in a s ln(l/ x) 
and in the classical field A% ~ 0(l/g). 

The mathematical structure of the CGC theory is rather complex and it will be only schematically 
described here. To that aim, it is convenient to switch to light-cone vector notations. Namely, for any 
4-vector such as x^, p^, A% etc. we shall define its light-cone (LC) components as 



1 

71 



(x° + x 3 ) 



1 

71 



X 



LC 



(26) 



In LC notations, the scalar product reads k • x = knX^ = k + x + k 



x 1 



k± • x±. 



To see the usefulness of these notations, consider a right-moving ultrarelativistic hadron, with 
(p^ o, 0, P) : this propagates at nearly the speed of light along the trajectory x 3 = t. In LC 
notations, the 4-momentum P^c — (v^^P, 0, 0, 0) has only a 'plus' component, while the trajectory 
reads simply x~ = 0. The same holds for any of the large-x partons which move quasi-collinearly 
with the hadron and serve as sources for the small-x gluons that we are interested in. In the semi- 
classical approximation, these small-x gluons are described as the solution to the Yang-Mills equations 
(the non-Abelian generalization of the Maxwell equations) having these 'fast' gluons as sources: 



DfF^(x) 



x ,x±) 



(27) 



In this equation, the l.h.s. features the covariant derivative D® b = d v — gf abc A u c and the field strength 
tensor = d v ' A% - d^A u a - gf abc A u b A% associated with the classical colour field, while the r.h.s. is 
the colour current of the 'fast' gluons: J a = p a , with p a (x~,x±) their colour charge density. The 
latter is localized in x~ near x~ = and is independent of time (hence of x + ), because these fast charges 
are 'frozen' by Lorentz time dilation. But the distribution of these charges in transverse space is random, 
since the fast gluons can be in any of the quantum configurations produced at the intermediate stages of 
the gluon evolution down to x. The proper way to describe this randomness is to give the probability 
to find a specific configuration p a {x~ , x±) of the colour charge density. This probability is a functional 
of p a (x~ known as the CGC weight function and denoted as Wy[p], with Y = ln(l/x). This 
functional is gauge-invariant, which in particular ensures that (p a (x~, x±)) = 0, as it should. 

To the accuracy of interest, all the observables relevant for the scattering off the small-x gluons 
are represented by gauge-invariant operators built with the classical field A%. If 0[A] is such an operator, 
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then its hadron expectation value is computed by averaging over all the configurations of p with the CGC 
weight function: 

(0[A]) Y = J Vp W Y [p] 0[A[p\] , (28) 



where A%[p] is the solution to Eq. (27). 



The expectation value ( [28] ) depends upon the rapidity Y = ln(l/x) via the corresponding de- 
pendence of the weight function Wy [p] • The latter is obtained by successively integrating the quantum 
gluon fluctuations in layers of x, down to the value of interest. One step in this evolution corresponds 
to the emission of a new gluon (with a probability 0(a s ) per unit rapidity) out of the preexisting ones. 
But unlike in the BFKL evolution, where gluons with different rapidities do not 'see' each other, in the 
context of the CGC evolution, the newly emitted gluon is allowed to interact with the strong colour field 



radiated by 'sources' (gluons and valence quarks) with higher values of x (see Fig. [14]). Accordingly, the 
change in the CGC weight function in one evolution step is non-linear in the background field Aa[p], 
and hence in the colour charge density p a . This procedure generates afunctional evolution equation for 
Wy [p] with the schematic form (see J8j[l0} for details) 

dWy[p] 



dY 



= H 



JIMWLK 



Wy[p], (29) 



where the JIMWLK Hamiltonian (from Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov, and Kovner 
fT8l[T9| ) i^jiMWLK is non-linear in p to all orders (thus encoding the rescattering effects in the emission 
vertex) but quadratic in the functional derivatives S/Sp (corresponding to the fact that there is only one 
new gluon emitted in each step in the evolution). In the dilute regime, where parametrically jp < 1, 
the non-linear effects are negligible, the JIMWLK Hamiltonian can be expanded to quadratic order in p, 
and then it describes the BFKL evolution. But for gp ~ 1, the non-linear effects encoded in //jimwlk 
prevent the emission of new gluons; this is gluon saturation. 

Eqs. ([27])— ([29]) are the central equations of the CGC effective theory. When completed with an 
initial condition at the rapidity Yq at which one starts the high-energy evolution, they fully specify the 
gluon distribution in the hadron wavefunction, including all its correlations. The initial condition Wy [p] 
is not determined by the effective theory itself, rather one must resort on some model. For a large nucleus 
(A » 1) and for Yq = 4 -=- 5 (corresponding to xo ~ 0.01), a reasonable initial condition is provided 
by the McLarren-Venugopalan (MV) model J7J, which assumes that the 'fast' colour sources are the 
N c x A valence quarks, which radiate independently from each other (since they are typically confined 
within different nucleons). The corresponding weight function is a Gaussian in p a . 

By taking a derivative w.r.t. Y in Eq. ( [28] ) and using Eq. ( |29| ) for Wy, one can deduce evolution 
equations for all the observables of interest. In general, these equations do not form a closed set; rather, 
they form an infinite hierarchy (originally derived by Balitsky |20|) which couples n-point functions with 
arbitrarily large values of n. In practice, this hierarchy can be truncated via mean field approximations 
(8][2T][22), leading to closed but non-linear equations, in particular the Balitsky-Kovchegov equation 
(24), that can be explicitly solved. It is also possible to numerically solve the functional JIMWLK 
equation ( [29] ), by first reformulating this as a stochastic process (a functional Langevin equation) 
which can be simulated on a lattice (25}|27j. 



In order to describe a scattering cross-section, the CGC effective theory developed so far must be 
combined with a factorization scheme. This will be described in the next subsection. 



3.5 Particle production from the CGC 

Let us start with some general remarks on factorization in scattering at high energies: this is a generic 
consequence of causality. For a hadron-hadron collision in the COM frame, the collision time At co n ~ 
1 / y/s is much shorter than the lifetime ([5]) of the partons participating in the collisions, which is propor- 
tional to the parton longitudinal momentum xP ~ y/s. Hence, these partons have been produced long 
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Fig. 15: Left: The space-time picture of a high-energy collisions, illustrating the factorization of the cross-section. 
Right: particle production in a proton-nucleus (or 'dilute-dense') collision. 



time before the collision, at a time where the two incoming hadrons were causally disconnected from 
each other (see Fig. [15] left). Accordingly, the respective parton distributions have evolved independently 
from each other and thus they are universal — i.e. independent of the scattering process that is used to 
probe them. This argument is purely kinematic and hence it remains true in the presence of QCD inter- 
actions leading to parton evolution or gluon saturation. However, the precise form of the factorization 
formula depends upon the kinematics and the structure of the process at hand and it is different when 
probing dense or dilute parts of the hadron wavefunction. 

(i) In the dilute regime, which corresponds to the situation where the transverse momenta p± of 
the produced partons are significantly larger than the saturation momenta in the two hadrons as evaluated 
at the relevant values of x, cf. Eq. fl3] ), the partonic subprocess involves merely a binary collision (cf. 
Fig. [6] left) : one parton in one projectile interacts with one parton in the other projectile, to produce the 
final state. Then, the cross-section depends only upon the parton densities (the 2-point correlations of 
the quark and gluon fields) in the incoming hadrons and the factorization formula takes a rather simple 
form: the hadronic cross-section is the convolution of two parton distribution functions (one for each 
hadron) times the cross-section for the partonic subprocess. 

Even in this case, one needs to distinguish between two types of factorizations, depending upon 
the kinematics of the final state: 



If the relevant values of x are not that small (say x > 0.01), then the parton evolution with 
decreasing x can be neglected and one can use the collinear factorization : the partons are assumed 
to move collinearly with the incoming hadrons (that is, one neglects their 'intrinsic' transverse 
momenta ~ Aqcd) and the parton distributions like xg(x, /j 2 ) depend only upon the longitudinal 
momentum fractions and upon the transverse resolution fj? (the 'factorization scale') of the hard, 
partonic, subprocess. The dependence upon /j 2 reflects the DGLAP evolution with increasing 
virtuality Q 2 . Eq. ( j_5 ) provides an example of collinear factorization. 

At smaller values of x, such that a s ln(l/x) > 1, the small-x evolution becomes important, lead- 
ing to an increase in the density of gluons and in their average transverse momentum. Because 
of that, the gluons cannot be considered as 'collinear' anymore : their distribution in transverse 
momenta must be explicitly taken into account. That is, one has to use the 'unintegrated' gluon 
distribution, whose evolution with 1/x is described by the BFKL equation, cf. Eq. (22). The 
corresponding factorization formula, known as hr-factorization, involves a convolution over the 
transverse momenta of the participating gluons. This is in fact a limiting case of the CGC factor- 
ization to be described shortly — namely its dilute limit, in which the saturation effects can be 
neglected. 
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(ii) The dense regime corresponds to collisions which probe saturation effects in at least one of 
the incoming hadrons. This happens when the transverse momenta of some of the produced particles 
are comparable with the saturation momentum at the relevant values of x. In such a case, the partons 
from one projectile scatter off a dense gluonic system (a CGC) in the other projectile, so they typically 
undergo multiple scattering. This is a non-linear effect similar to saturation: each additional scattering 
represents a correction of order a s n to the cross-section, which for n ~ l/a s is an effect of order 
one. (As usual, n = k±) denotes the gluon occupation number in the dense projectile and is of 
0(l/a s ) when k± < Q s (x).) So, when a parton scatters off a CGC, the multiple scattering series 
must be resummed to all orders. This resummation involves arbitrarily many insertions of the strong 
colour field A% which represents the CGC, which implies that the associated cross-section is sensitive to 
multi-gluon correlations (n-point functions of the field A% with n > 2). Clearly, the multiple scattering 
cannot be encoded in the (collinear of kr) factorization schemes alluded to above, which involve only 
the respective 2-point functions — the parton distributions. 

Fortunately, there is an important simplification which occurs at high energy and which permits 
to compute multiple scattering to all orders: an energetic parton is not significantly deflected by the 
scattering and thus can be assumed to preserve a straight line trajectory throughout the collision. This 
is known as the eikonal approximation. To explain it in a simple, but phenomenologically relevant, 
setting, consider first proton-nucleus (p+A) collisions, cf. Fig. [15] right. This is an example of a dense- 
dilute scattering, that is, a collision in which a projectile which is relatively dilute (the 'proton') and can 
therefore be treated in the collinear or the kr factorization, scatters off a dense target (the 'nucleus'). 
Then the partons from the dilute projectile undergo multiple scattering off the strong colour field of 
the dense target. For the kinematical conditions at RHIC or the LHC, this 'dense-dilute' scenario is 
optimally realized in the case of particle production at forward rapidities, that is, in the fragmentation 
region of the proton (deuteron at RHIC). 

To be specific, consider the production of a light quark with rapidity y > and semi-hard trans- 
verse momentum p±. The production mechanism is as follows: a quark from the proton, with relatively 
large longitudinal momentum fraction x\ = (p±/\fs)e y and negligible transverse momentum, scatters 
off the gluons with X2 — (p±/\fs)e~ y and &2_L ~ Qs(%2) from the nucleus — which form a dense sys- 
tem because X2 <C 1 and A » 1 — and thus accumulates a final transverse momentum p± ~ Q s {x2)> 
Within the CGC effective theory, the nucleus in a given scattering event is described as a classical colour 
field A%, off which the quark scatters with the ^-matrix 

S aP = (P\ Texp [i J dVintG/)}|«> with £ [nt (y) = f a {y)A^y) and f a {y) = q^H^ . 

(30) 

Here Ant (y) is the Lagrangian density for the interaction between the colour current ja of the quark and 
the colour field A^ of the target, and a (/3) represents the ensemble of the quantum numbers character- 
izing the state of the quark prior to (after) the collision. The quark deflection angle reads 9 ~ p±/Ei 
with Ei = xiyfs/2 (the quark energy). For the kinematics of interest we have p± <C Ei, so this an- 
gle is small, 9 <C 1, as anticipated. Hence, one can assume that the quark keeps a fixed transverse 
coordinate x± while crossing the nucleus. This is the eikonal approximation. In this approximation, 
a = (x±,i) and (3 = (x±,j), where i and j are colour indices in the fundamental representation of 
SU(iV c ) which indicate the quark colour states before and respectively after the scattering. Also, assum- 
ing the quark to be a left-mover (and hence the nucleus to be a right mover, as in Eq. ([27])), one can write 
j%{y) ~ 5v-gt a 5(y+)5( 2 \y ± - x±). Then Eq. ^ reduces to 

{x ± J\S\x ± ,i) ~ Vijix^) with V(x ± ) =Texp{i 5 J dx~A+(x-, x ± )t a } , (31) 

where the 'path-ordering' symbol T denotes the ordering of the colour matrices A+(x~, x±)t a in the 
exponent, from right to left, in increasing order of their x~ arguments. The integration runs formally 
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over all the values of x~, but in reality this is restricted to the longitudinal extent of the nucleus, which 
is localized near x~ = because of Lorentz contractioij^] The path-ordered exponential V is a colour 
matrix in the fundamental representation, also known as a Wilson line. It shows that the only effect of the 
scattering in the high energy limit is to 'rotate ' the colour state of the quark while the latter is crossing the 
nucleus. If instead of a quark, one would consider the scattering of a gluon, the corresponding S-matrix 
would be again a Wilson line, but in the adjoint representation (t a — » T a ). When the target field is weak, 



gA + <C 1, one can expand the exponential in Eq. (31 ) in powers of gA + , thus generating the multiple 



scattering series. But when gA + ~ 1, as is the case for a target where the gluons are at saturation, such 
an expansion becomes useless, since all the terms count on the same order. In such a case, one has to 
work with the all-order result, as compactly encoded in the Wilson line. 

The above considerations also show that multiple scattering at high energy is most conveniently 
treated in the transverse coordinate representation: the successive collisions modify the transverse mo- 
mentum of the partonic projectile, but do not significantly alter its transverse coordinate (or 'impact 
parameter'). But the interesting observable is the cross-section for producing a quark (or gluon) with a 
given transverse momentum p± and rapidity y. This is obtained by multiplying the amplitude Vij(x±) 
with the complex conjugate amplitude Vj i (y±) for a quark at a different impact parameter y±, and then 
taking the Fourier transform x± — y± p±. This yields (for the forward kinematics of interest here) 



~ xrf^xupl) J d 2 r ± e~ ir ^ ±- (trV(x 1 _)V^y 1 _)) Y . (32) 



dy d 2 p A 

The quark distribution xif q {x\ 1 p 2 A ) gives the probability to find a quark 'collinear' with the proton, with 
longitudinal momentum fraction x\, on the resolution scale p\ set by the partonic scattering. Within the 
integral, we have defined r± = x± — y±. Furthermore, the colour trace has been generated by summing 
over the final colour indices (J2j) anc ^ averaging over the initial ones ((1/N C ) ^), and the brackets 
denote the average over the target wavefunction, evolved up to the rapidity Y = In (1/^2). In the CGC 
formalism, this target average is computed according to Eq. ([28]), that is, as an average over the colour 



charge density of the 'fast' sources which are responsible for the target field {x , x±) via Eq. (27): 



S Y (x ± ,y±) = ^- (tvV(x ± )V\y ± )) Y = f VpW Y [p] ^;tr(V(s L )W(y_ L )) . (33) 

This 2-point function of the Wilson lines is recognized as the S'-matrix for the scattering between a 
colour dipole (a quark-antiquark pair in an overall colour singlet state) and the CGC. The non-linear 
effects enter this S'-matrix at two levels: (i) via multiple scattering for the quark and the antiquark, 
as described by the respective Wilson lines, and (ii) via gluon saturation in the target wavefunction, 
as encoded in the CGC weight function. In this context, the target saturation momentum Q S (A,Y) 
also plays the role of the unitarization scale for the colour dipole: the dipole scattering becomes strong 
(meaning that |Sy(rj_)| <C 1) when the dipole size r± is of order l/Q s or larger. Indeed, so long as 
the dipole is relatively small, such that r± <C it predominantly scatters off the gluon modes with 

k± ~ l/r_L ^> Q s , which are dilute. The corresponding target field is weak (gA + <C 1), the Wilson 
lines are close to one, and so is the dipole S-matrix. Namely, for r±Q s <^ 1 one finds 1 — Sy(r^) oc 
(r[(Js(y)) 7s with 7 S ~ 0.63. On the other hand, a large dipole with r± > 1/Q S (Y) probes the high- 
density gluon modes with k± < Q S (Y) \ the associated colour fields are strong, gA + ~ 1, so the Wilson 
lines are rapidly oscillating and their product averages out to a very small value: \Sy(t±)\ <C 1 when 
r ± Q s > 1. 

The dipole S-matrix provides a convenient framework to study high energy evolution and sat- 
uration since, in the limit where the number of colours is large 7V C » 1, it obeys a relatively simple 
equation for the evolution with Y — a non-linear generalization of the BFKL equation known as the 



More precisely, the small-x gluons which participate in the scattering are delocalized within a distance Ax ~ 1/ (X2-P) 



around x = 0, as explained after Eq. ( 18 1. 
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Fig. 16: The results of the JIMWLK evolution for the generalized gluon distribution in Eq. (35), parameterized 
as <&(Y,k\_) = (l/a s )C(Y : k±). Left: as a function ofk±L, with L the transverse size of the system. Right: as 
a function of the scaling variable k±/Q s (Y), to exhibit geometric scaling (the curves corresponding to different 
values ofY fall approximately on top of each other). From Ref. [26]. 



Balitsky-Kovchegov (BK) equation (24). Originally deduced within Mueller's 'dipole picture' (9][28) 
(an insightful reformulation of the BFKL evolution valid at large N c ), the BK equation also emerges 
from the large- N c limit of the Balit sky- JIMWLK hierarchy. Remarkably, the BK equation is presently 
known to next-to-leading order (NLO) accuracy [29}|3T| . This is important in view of phenomenologi- 
cal studies of both deep inelastic scattering — at high energy, the DIS cross-section ([9]) can be related to 
the dipole S-matrix (32) and the NLO corrections are essential in order to achieve a good description of 
the HERA data at small x (33] - [37) — and of particle production in dense-dilute scattering, which is the 
topics of interest for us here. 

Specifically, Eqs. ([32])— ([33]) express the CGC factorization for inclusive quark production in dense- 



dilute scattering fTTj[38]-[40l|. Eq. (32) is usually written as 



c ^x x f q {xuA)*A<y,p L ,b L ), (34) 

where the first factor a s jp\ is the cross-section for the elementary q + g — >> q scattering (the would-be 
partonic subprocess in the single scattering limit) while 

2 r i 

$A(y,P±,&i) = — / d 2 r ± e~ ir ^ — (tiV(x ± )V\y ± )) , (35) 
a s J N c \ Iy 

plays the role of a generalized unintegrated gluon distribution (here, for the nucleus) at impact parameter 



b± = (x± + y±)/2. When p± » Q s (^4, Y), Eq. (35 ) can be evaluated in the single scattering approx 



imation, as obtained by expanding the product of Wilson lines to quadratic order in A + . In that (dilute) 
regime, reduces indeed to the usual unintegrated gluon distribution — the one which enters the fcy- 
factorization and obeys the BFKL equation. But for lower momenta p± < Q S (A, Y), the non-linear 
effects become essential, as already discussed in relation with the dipole scattering. 



Fig. [16] shows the CGC prediction for the generalized gluon distribution ( [35] ), as obtained via 
the numerical resolution of the JIMWLK equation with initial conditions at Y = of the McLerran- 
Venugopalan type (and with a running coupling) (26). The two plots exhibit a pronounced peak at a 
special value of the transverse momentum p± which increases with Y : this special value is, of course, 
the saturation momentum Q S (A, Y). In fact, it is quite easy to understand these plots in the light of 
our previous discussion of the dipole scattering. The Fourier transform in Eq. ( [35] ) is controlled by the 
competition between the complex phase e~ ir± p± and the dipole S-matrix Sy(r±). For large momenta 

27 



P± ^> Qs, the complex exponential limits the integration to small dipole sizes r± < l/p± <C for 
which the scattering is weak: 1 — Sy (r±) ~ ( r j_QsQO) 7s . Then the Fourier transform yields 



$a(Y, P± ) ~ - f g ' ( f r) V 5 for p ± » Q,(A,y). 
a» V P± ) 



(36) 



pi 

The difference 1 — 7 S ~ 0.37 is an anomalous dimension introduced by the high energy evolution. 



(Without this evolution, one would have the bremsstrahlung spectrum oc cf. Eq. (22).) For 

lower momenta p± < Q s , the integral over r± in Eq. ( |35] ) is limited by the dipole S'-matrix to values 
r± < 1/Qs < Vpl- (Recall that S Y (r±) < 1 when rj_ > 1/Q S .) Then 



Z P ~\ Z 

<S>a(Y, P± ) ~ ^ / d 2 r ± 9(1/Q 8 - r ± ) ~ — P j- — for p ± < Q S (A, Y) . (37) 
Eqs. ( [36] ) and ( [37] ) explain the pronounced peak in the gluon distribution at p± ~ Q S (Y), as visible 



in Fig. 16 This in turn implies that the cross-section for particle production is dominated by semi- 
hard gluons with p± ~ Q S (X) and hence can be computed in perturbation theory. Another important 
consequence of saturation, which is visible too in Eqs. ( [36] ) and ( [37] ) (and is numerically tested in the 



right panel of Fig. [16]), is geometric scaling pT| - [44[ : the unintegrated gluon distribution depends upon 
the two kinematical variables p± and Y only via the dimensionless ratio pj_/Q s (Y). This scaling has 
important consequences for the phenomenology, to be discussed in the next section. 



Note also that, as a result of the non-linear physics, the generalized 'gluon distribution' ( [35] ) is a 
process-dependent quantity : it depends not only upon the gluon density in the target, but also upon the 
nature of the partonic subcollision. For instance, if instead of the quark production, one would consider 
the production of a gluon, the 'fundamental' Wilson lines in Eq. ( [35] ) would be replaced by 'adjoint' ones 
| [38] - [40| . Also, if one considers a more complicated final state — say, the production of a pair of partons 
— then the analog of Eq. ( [35] ) will involve a pair of Wilson lines for each of the partons partaking in the 
collision (one such a line in the direct amplitude and another one in the complex conjugate amplitude). 



More examples in that sense can be found in Refs. |45]-[50|. See also (5TJ for a recent, comprehensive, 



discussion of the relation between the CGC factorization and the /cT-factorization. These considerations 
show that, in the presence of non-linear effects, the notion of 'parton distribution' ceases to be useful: 
the observables involve higher n-point functions of the gluon fields (generally, via the product of Wilson 
lines) which moreover couple with each other under the non-linear, JIMWLK, evolution. The complete 
information about the gluon correlations and their evolution with Y (to leading-logarithmic accuracy at 
least) is encoded in the CGC weight function. 

We now turn to nucleus-nucleus (A + A) collisions, which for the typical kinematical condi- 
tions at RHIC and the LHC represents an example of dense-dense scattering : the wavefunctions of 
both nuclei develop saturation effects which influence the production of particles with semi-hard trans- 
verse momenta. Hence, in order to compute the bulk of particle production, one must take into account 
the many-body correlations associated with gluon saturation in both nuclei, together with the multiple 
scattering between these two saturated gluon distributions. Once again, this complex problem can be ad- 
dressed within the CGC formalism and in the eikonal approximation. Not surprisingly, the treatment of 
the two nuclear projectiles is now symmetric, they are both described as colour glass condensates, with 
weight functions Wy x [p\] and, respectively, Wy 2 [p2]- And their collision is described as the scattering 
between two classical distributions of colour charges moving against each other. 

To be more specific, consider inclusive gluon production in the COM frame, where the nuclei have 
rapidities ±?/beam- If the produced gluon has rapidity y, it will probe the evolutions of the two nuclear 
wavefunctions up to rapidities Y\ = ^beam — V and Y2 = 2/beam + V, respectively. This motivates the 
following CGC factorization for the spectrum of the produced gluons (52) 

(38) 

class 
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Fig. 17: Left panel: a cartoon of the 'dense-dense' collision between two heavy nuclei, which illustrates the 
complexity of the process. Non-linear effects enter via gluon saturation (in both nuclei) and multiple scattering. 
Right panel: the gluon spectrum generated by this collision within the CGC formalism (from Ref. f26]/). 



where the last factor inside the integrand, (d7V/d?/d 2 p^)| c i ass , represents the spectrum produced in the 
scattering between two given configurations of classical colour sources (the 'fast partons') — one for 



each nucleus. More precisely, as discussed in relation with Eq. ( |27| ), the right-moving nucleus is de- 
scribed in a given event as a colour current having only a 'plus' component: Jf' a = 5^ p\ , with the 
charge density p\ localized near x~ = (due to Lorentz contraction) and independent of x + (by Lorentz 
time dilation). Similarly, the left-moving nucleus is represented by a colour current with only a 'minus' 
component: J^' a = p§ ? with p\ localized near x + = and independent of x~. At a classical level, 
the 'scattering' between these two currents is described by the solution A% to the Yang-Mills equation 
including both types of sources (compare to Eq. ( |27| )) : 

DfF^(x) = 6" + pUx-,x ± ) + 6*-(%(x + ,x 1 _). (39) 

This equation describes multiple scattering because it is non-linear: the collision begins at x + = x~ = 
(i.e. t = z = 0) and for positive values of x + and x~, the solution A% is non-linear to all orders 
in both pi and p2. This solution cannot be computed analytically, but numerical solutions are by now 
available J26{[53] - [55| . The cross-section (d7V/dyd 2 p^)| c i ass for particle production is obtained via the 
Fourier transform of this classical solution, that is, by projecting the field A% onto modes with transverse 



momentum p±. Finally, the average over the CGC weight functions of the two nuclei, cf. Eq. ( [38] ), is 
numerically performed. It is this last procedure which introduces the dependence of the cross-section 
upon the rapidity y, via the corresponding dependence of the two weight functions. Note that, in line with 
the general philosophy of the CGC formalism, the only quantum effects to be included in the calculation 
are those associated with the high-energy evolution of the projectile wavefunctions prior to scattering, 
which are enhanced by the large logarithms Y{ = ln(l/x^), with i = 1, 2. The final outcome of this 
calculation is the gluon spectrum displayed in Fig. [17] (right panel) (26). This is very similar to the 
'unintegrated' gluon distribution in any of the incoming nuclei (compare to Fig. [16]), in particular, it is 
peaked at a value of p± of the order of the saturation momentum and which increases with y. Some 



further consequences of the solution to Eq. ( p9| ) will be discussed in Section [377. 

Note finally an important difference between p+A (dilute-dense) and A+A (dense-dense) colli- 
sions: in the former, the particles produced by the collision do not interact with each other, but merely 
evolve via fragmentation and hadronisation towards the final hadrons observed in the detectors; by con- 
trast, in A+A collisions the partonic medium created in the early stages of the collisions is very dense, 
so these partons keep interacting with each other — one then speaks about final state interactions, as 
opposed to the initial state ones, which were associated with high density effects like saturation in the 
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incoming wavef unctions. These 'final state' interactions redistribute the partons in energy and momen- 
tum, which makes it difficult to verify the early-stages spectrum, as given by the aforementioned CGC 
calculations, against the measured hadron yield. Yet, the CGC formalism has the virtue to provide the 
initial conditions for the subsequent dynamics, at a time r ~ 1/Q S ~ 0.2 fm/c. So, its predictions can 
be at least indirectly tested, via calculations of the final state effects which include initial conditions of 
the CGC type. We shall return to such issues later on. 



3.6 Some experimental signatures of the CGC in HIC 

Let us consider now some phenomenological applications of the previous results for 'dilute-dense' scat- 
tering |T^[TT}[T4j[T5j (and Refs. therein). We start with the R p a ratio, defined as the ratio between 
particle production in p+A collisions and that in p+p collisions for the same kinematics; schematically, 



R p a(v,p±) = 



NcoU 







d 2 p_\_dr] 


pA 


dN h 




d 2 p_\_dr] 


pp 



(40) 



where the subscript h denotes the hadron species and N co u is the number of binary proton-nucleon 
collisions in the p+A scattering at a given impact parameter, as computed under the assumption that 
the various nucleons inside the nucleus scatter independently from each other. For relatively central 



collisions (i.e. small impact parameters; see Fig. 23), one has N co u ~ A 1 / 3 . The normalization in 



Eq. (40 ) is such that R p a would be equal to one if the p+A collision was a superposition of A incoherent 
p+p collisions. Conversely, any deviation in R p a from unity is an indication of coherence (high-density) 
effects in the nuclear wavefunction. Such a deviation is clearly seen in the respective RHIC data at 
forward rapidities (77 > 0). More precisely, RHIC performed deuteron-gold (d+Au) collision^] at = 
200 GeV per nucleon pair and measured the ratio i?d+Au f° r semi-hard momenta p± = 1 5 GeV and 



for rapidities 77 = 4 in the deuteron fragmentation region (56}[57 1 



The results of the corresponding analysis by BRAHMS are shown in Fig. 18 For p± > 2 GeV, 
they show an enhancement (i?d+Au > 1) f° r 77 = 0, known as the 'Cronin peak', which however 
disappears when increasing 77, leading to suppression (i?d+Au < 1) at 77 > 1. E.g. for 77 = 3.2, one finds 
^d+Au — 0.6-1-0.8 for p± = 2-^4 GeV. This behaviour can be understood in terms of saturation in the 
nuclear gluon distribution and its evolution with Y (58]-[6TJ. For central rapidities 77 ~ andp^ = 2 GeV, 



one probes X2 ~ 10 -2 , which is large enough for the high-energy evolution to be negligible. Then the 
gluon density in the nucleus is large juste because of the many (3A) valence quarks acting as sources 
for gluon radiation. An incoming parton from the deuteron scatters off this dense gluonic system and 
thus acquires an additional transverse momentum of order Q S (A, X2) ~ 1 GeV. This yields a shift in the 
spectrum of the produced particles towards higher values of p±, leading to the Cronin peak. For forward 
rapidities, say 77 = 3.2 and p± = 2 GeV, one has X2 ~ 10 -4 and then the high-energy evolution with 



Y = In (1/^2) is important. In that case, one can use the factorization ( [34] ) for both d+Au and for the 
p+p collision which serves as a benchmark, leading to 

^ 1 $ A (Y, P± ) 
Rd+A —A^%(Y^j- (41) 

For this kinematics, the saturation effects are important in the gold nucleus (since the nuclear saturation 
scale Q S (A, Y) ~ 2 GeV is comparable with p±) and lead to a slow down of the evolution with Y. 
On the other hand, saturation is still negligible for the proton, so the corresponding unintegrated gluon 
distribution <f> p (Y,p±) rises rapidly with Y, according to the BFKL evolution. Hence, when increasing 
77 (and thus Y), the denominator in Eq. ( |4T| ) rises much faster than the numerator there, leading to a 
decrease in the ratio. This is precisely the trend seen in the data. 



5 For d+Au collisions the number of binary collisions in Eq. ( |40| ) should be evaluated as N co ii — 2 A 1 7/3 since the projectile 
deuteron involves 2 nucleons, i.e. twice as much as the proton. 
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Fig. 18: BRAHMS results for d+Au collisions at 200 GeV/nucleon f56| / ; the ratio RdAu for charged hadrons as a 
function ofp± for central and forward pseudorapidities. 



In particular, for sufficiently high Y, such that both the nucleus and the proton in the ratio in 
Eq. ( |4T] ) are at saturation, and for transverse momenta comparable or sli ghtly larger to the nuclear satu- 



R d+ An * ~JT7s ( ^P^V ) * A- 1 ^ , (42) 



ration momentum (p± > Q S {A, Y) ~ A 1 / 3 Q s (p, Y)), one can use Eq. (36) for the unintegrated gluon 
distributions in both targets. This yields a particularly simple result, 

1 ' Ql(A,Y) ^' 
AV3 \Q*(p,Y) 

which is independent of both the transverse momentum (within a limited range above the nuclear satu- 
ration momentum) and the rapidity. It would be interesting to test this prediction in p+A collisions at the 
LHC, where the coverage in both Y and p± will be larger than at RHIC. 

Another important consequence of saturation, which is manifest in Eqs. ([36])— ([37]) and has been 



implicitly used in deriving Eq. ( |42| ), is geometric scaling in the gluon distribution. Via Eq. ( |34| ), this 
suggests a similar scaling in the spectrum of the produced particles. This prediction has been tested 
against the LHC data for p+p collisions (62), with the results shown in Fig. [19] There one can see the 
ratio 

(djv/dyd„)| 

= pssic <43) 

between the measured spectra for single-inclusive charged hadron production at two COM energies, 
Ei = y/siandE2 = y/s2, and for midrapidities: \rj\ < 2.4. More precisely, one displays two such ratios, 
as obtained by combining the LHC data for three different energies: y/s = 0.9, 2.36, and 7 GeV. (Note 
that these energies are high enough for the saturation effects to be important in the proton wavefunction 
for 7] ~ and semi-hard transverse momenta.) If the spectrum (dN /d 2 p±dr]) scales as a function 
of r = p±/Q s (Y), then the ratio ( [43] ) should be equal to one when plotted as a function of r. This 
expectation is indeed met by the data, as shown in the r.h.s of Fig. [19] At this point, it is worth noting 



that geometric scaling has been first observed in the HERA data for DIS (4TJ and that the experimental 
search for this remarkable behaviour has been inspired by the theoretical ideas about gluon saturation. 

Another remarkable regularity in the data which is naturally explained by the CGC is the limiting 
fragmentation : when the charged particle rapidity distribution dN/ dr] is plotted as a function of the 
variable r/ = 77 — ybeam — the rapidity difference between the produced particle and the dilute projectile 
(say, the deuteron in the case of d+Au collisions) — , then the distribution turns out to be independent of 



the collision energy over a wide range around 7/ = 0, whose extent is increasing with (see Fig. 20) 



Note that rj f ~ corresponds to forward rapidities (77 > 0) according to our previous terminology, to 



which Eq. (34) applies. Moreover, for such rapidities, even a nucleus-nucleus collision may be viewed 
as 'dilute-dense', in the sense that one of the nuclei is probed at large x\, where it looks dilute. 

To understand limiting fragmentation on the basis of Eq. (|34|), notice that (i) 77 — ybeam — l n ^i 



(the rapidity of the produced particle is roughly equal to that of the fast parton which initiated the scatter- 
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Fig. 19: The ratio Re 1 /e 2 between particle spectra at energies E\ and E 2 , as measured for three energies at the 
LHC, is plotted as a function ofp± (left panel) and of the scaling variable r = p±/Q s (Y) (right panel). From 
Ref. $62$. 




Fig. 20: Pseudorapidity density distributions of charged particles emitted in d+Au collisions (left panel) and 
Au+Au collisions (right panel, separately for central and peripheral impact parameters), at various energies, as a 
function of the rj' = r\ — 2/beam variable. In the right panel, the impact parameter is represented by the centrality 
bins: 0-6% (i.e. the 6% most central data) and respectively 35-40%. From Ref. %64§ . 



ing), (ii) the 'multiplicity' dN/dr] * s obtained by integrating the spectrum ( |34| ) over all values of p±, 
(iii) this integral is dominated by p± ~ Q S (Y) (since the gluon distribution is strongly peaked at Q s , 
cf. Fig. 16), and (iv) the result of the integration is very weakly dependent upon Y = ln(l/x2), by ge- 
ometric scaling (indeed, <&a(Y,p±) is roughly a function of r = p±/Q s (Y), which is evaluated at r = 1 
when performing the integral). These arguments imply that dN/dr] depends upon x\ ~ exp(?7 — ybeam) 
but is approximately independent of X2 (and hence of the total collision energy) within the range of 'for- 
ward rapidities'. This is precisely the property called limiting fragmentation, as visible in Fig. |20| for 
both d+Au and Au+Au collisions at RHIC (PHOBOS) (63j[64). 

As a final application for 'dilute-dense' scattering, we shall consider the production of a pair of 
hadrons at forward rapidities. The kinematics has been already explained in relation with Eq. ( [T3] ) : 
when both r\ a and 775 are positive and large, the produced hadrons explore small values X2 <C 1 in the 
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Fig. 21: Di-hadron azimuthal correlations at forward rapidities (rj a , rjb ~ 3) and semi-hard transverse momenta 
as measured in p+p (left panel) and d+Au (right panel) collisions at RHIC, for V200 GeV per nucleon pair. 



target wavefunction and hence they can experience high density effects — multiple scattering and gluon 
saturation. As before, these effects are important when the transverse momenta p a ± and of the 
two hadrons are comparable to the target saturation momentum. The respective cross-section admits a 
factorization similar to Eq. ( [32] ), where however the target expectation value now involves the trace of the 
product of four Wilson lines: two for the produced partons in the direct amplitude, and two for the same 
partons in the complex conjugate amplitude. Hence, the generalized 'gluon distribution' of Eq. ( [35] ) gets 
now replaced by a 4-point function, known as a colour quadrupole [22, 47] - [50| . 

A convenient way to study high-density effects in this setup is to measure di-hadron correlations 
in the azimuthal angle, which is the angle indicating the direction of propagation in the transverse plane. 
If the medium effects are negligible, two relatively hard particles (pj_ » Aqcd) are produced back-to- 
back (p a A_ + Pb± — 0). So, if the trigger detects one of these particles together with its fragmentation 
products (a 'jet'), then by measuring the particle distribution in the same event one should find another 
'jet' at a relative angle A<3> 2^ tt. On the other hand, if the target looks dense on the transverse resolution 
scale of the produced particles, then the ^-distribution gets broaden via multiple scattering and the peak 
corresponding to the 'away jet' at A<3> 2^ tt gets smeared, or it even disappears. 

This scenario has been indeed confirmed by the RHIC data as measured by STAR [65 ,66]. Fig. [21 



shows the experimental results for di-hadron production at forward rapidities (j] a ^ ~ 3) and semi-hard 
transverse momenta in both p+p and d+Au collisions. As already discussed in relation with Eq. ( |40| ), for 
this kinematics one expects the saturation effects to be negligible for a proton target, but important for a 
gold nucleus. And indeed, the data for azimuthal correlations in Fig. [21] show a pronounced 'away' peak 
in the p+p collisions (the left panel), but a strongly suppressed one in the d+Au ones (the right panel). 
Such a suppression was actually predicted by the CGC effective theory (49}[67j and its experimental 
observation at RHIC is one the most compelling evidences in favour of gluon saturation available so far. 

Turning to A+A (or 'dense-dense') collisions, we recall that, in that case, the CGC framework 
provides the initial conditions for the subsequent evolution of the liberated partonic matter, but not also 
the spectrum of the produced hadrons. Still, predictions can be more for more inclusive quantities like 
the total multiplicity, which are expected to be less affected by 'final state' interactions. In particular, 
Fig. [22] exhibits the centrality dependence of the multiplicity of the charged particles as measured at the 
LHC in Pb+Pb collisions at y/s = 2.76 TeV/nucleon pair, together with the respective predictions of 
various theoretical models. The centrality of a collision refers to the relative impact parameter b± of 
the two projectiles in the transverse plane (see Fig. [23]). For A+A collisions, this is often parameterized 
in terms of the 'number of participants' — the number of incoming nucleons from the two nuclei in 
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Fig. 22: The LHC data (ALICE ^6E\j ) for the midrapidity (rj = 0) charged particle multiplicity in Pb+Pb collisions 
at y/s = 2.76 TeV/nucleon pair, normalized by the number of participants N part and plotted as a function of N part 
— hence of the centrality of the collisions, which increases (meaning that the impact parameter b± decreases) with 
increasing N part . The predictions of six theoretical calculations are also shown. The three lowest ones, marked 
as 'Armesto et al\ 'Kharzeev et aV , and Albacete et aY , are inspired by the CGC effective theory and explicitly 
include gluon saturation. The three other models include some of the effects of saturation, in the form of an 
energy-dependent 'infrared' cutoff at low transverse momenta, p™ in ~ E x , which effectively plays the role ofQ s . 



the region where the nuclei overlap with each other (the 'interaction region'). Clearly, central collisions 
~ 0) involve more participants than the peripheral (large b±) ones. Although one cannot compute 
d7V c h/d?7 fully from first principles (as this also requires some information about the distribution of 
nucleons within the nuclear disk), one can easily estimate the dependence of this quantity upon iV par t 



within the framework of the CGC effective theory. As discussed in relation with Fig. [17} the partons 
produced in the early stages of a A+A collision are typically gluons with transverse momenta p± ~ Q s , 
which have been liberated by the collision. So the multiplicity diV^/dr/ near 77 = is proportional to 
the number of such gluons which were present in the initial nuclei, within their region of overlapping: 



dN ch 



xg A (x,Q 2 s ) _ 1 Q ~ fA v AT _ X .Q 2 8 (A,E) 



o, S * AV oc ——SQj(A,E) ~ AUt£ A ln ; v 2 ' 7 . (44) 

drj 77=0 nR 2 A a s (Q 2 s ) Aqcd 



Here, S is the transverse area of the interaction region and we have used Eq. ( |23| ) for the (nuclear) sat- 
uration momentum together with the fact that S Q 2 (A, E) is proportional to iVpart an d it grows with 
the COM energy E = like E x (cf. Eq. (25) and the discussion of Fig. [l3]). Hence, the ratio 



(dAf c h/d77)/Af part is expected to be only weakly dependent upon Af part , via the corresponding depen- 
dence of the running coupling: l/a s {Q 2 ) ~ ln(Q^ / ^Aq CD ) ~ In A^ part . This is in good agreement with 
the data, as shown in Fig. |22j (For the most refined calculation to date, whose results are indicated in 
Fig.|22]by the curve denoted as 'Albacete et al', see Ref. (69).) 

3.7 The Glasma 

By the uncertainty principle, it takes a time At ~ l/Q s to liberate a particle with transverse momentum 
P± ~ Qs at midrapidities (77 ~ 0). So, by that time, the small-x gluons which were originally confined 
within the wavefunctions of the incoming nuclei, are released by the collision. What is the subsequent 
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Fig. 23: A sketch of the geometry of a heavy ion collision. The collision (or longitudinal) axis is denoted as 
z, while x and y are the transverse coordinates. The interaction region (the almond-shape region where the two 
nuclei overlap which each other) is singled out. N part is the number of nucleons in this region. The interaction 
region is horizontally cut by the reaction plane (x,z). 



evolution of these gluons ? To answer this question one needs a better understanding of their configura- 
tion at the time of emission. This can be inferred from the solution A% to the Yang-Mills equation ([39]), 
or, more precisely of the associated chromo-electric and magnetic fields, E a and B a . 

Prior to the collision (t < 0), these fields describe the two 'colour glass condensates' of the 
incoming nuclei, which by causality are independent of each other. Due to the high-energy kinematics, 



the CGC fields turn out to be quite simple (see Fig. [24]) : for each nucleus, the respective vectors E a 
and B a have only transverse components, E l a and B l a with i = 1, 2, meaning that they are orthogonal 
to the collision axis x 3 . Besides, they are also orthogonal to each other, E a • B a = (for each of the 
N% — 1 = 8 values of the colour index a), and they have equal magnitudes: 

E a JL B a JL x 3 , \E a \ = \B a \ (prior to the collision) . (45) 

These initial electric and magnetic fields are localized near x~ = for the left-moving nucleus and, 
respectively, x + = for the right-moving one, like the respective colour charges. In a given event, 
their values and orientations can randomly vary from one point to the other in the transverse plane. 
But on the average, the fields at different points are correlated due to 'memory' effects in the high 
energy evolution, in particular, due to saturation. The correlations, which are encoded in the respective 
CGC weight functions, are typically restricted to a saturation disk, i.e. to transverse areas with radius 
~ 1/Qs : domains separated by transverse distances Ax± ^> 1/Q S evolve independently from each 
other, since saturation prohibits the emission of gluons with momenta k± <C Q s . Within a saturation 
disk, gluons arrange themselves in such a way to shield their colour charges and thus minimize their 
mutual repulsion; accordingly, a saturation disk has zero overall colour charge (see the right figure in 
Fig. [3]). Also, gluons can be correlated with each other in rapidity, due to the fact that they have common 
ancestors, i.e. they belong to the same parton cascade. Such correlations extend over a rapidity interval 
Ay rsj l/a s , since this is the typical value of Y which is required to build parton cascades according 
to the BFKL evolution. These correlations, which were built in the initial wavefunctions via the high- 
energy evolution, get transmitted to the gluons liberated by the collision and thus have consequences for 
the distribution of the particles in the final state. 

In view of the above, a A+A collision can be viewed as the scattering between two sheets of 



coloured glass, as illustrated in Fig. 24 Incidentally, a similar structure for the incoming fields — electric 



and magnetic fields which are orthogonal to the beam axis and to each other, and which are localized 
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Fig. 24: The chromo-electric and chromo-magnetic field configurations in the two nuclei prior to the collision. 



near the respective light-cone and frozen by Lorentz time dilation — would also hold if the nuclei were 
made with electric (rather than colour) charges. In both QED and QCD, such field configurations — 
known as the Weizsacker-Williams fields — represent the boosted version of the Coulomb fields created 
by the ensemble of (electric or colour) charges in their rest frame. However, the non-Abelian structure 



of QCD is essential for having a collision: the Abelian version of Eq. (39) would be linear and it would 



not describe a scattering process. (The total electromagnetic field at t > would be simply the sum of 



the individual fields of the two nuclei.) The non-linear effects encoded in the Yang-Mills equation p9y 
describe the scattering between the small-x gluons in the two CGC's. The corresponding solution at 
t > (more precisely, in the forward light cone at x + > and x~ > 0, which is the space-time region 
causally connected to the collision) represents the gluonic matter produced by the collision. 

This solution exhibits a very interesting structure: in addition to the transverse fields on the two 
sheets, which after the collision are separating from each other, there are also longitudinal, electric and 
magnetic fields, E 3 and B 3 , which extend along the collision axis. The latter give rise to colour flux tubes 



(or 'strings') with the endpoints on the two sheets and a typical transverse radius l/Q s (see Fig. 25 left). 
Right after the collision (t = + ), these fields are quite strong, E 3 ~ B 3 ~ 1/g, since they carry most 
of the energy of the original CGC fields. At such early times, the gluonic matter is still in a high-density, 
coherent state, for which a description in terms of classical fields is better suited than one in terms of 
particles. But with increasing time, the system expands, its density decreases, and so does the strength of 
the fields. After a time t ~ 1/Q S , the magnitudes of all the fields (transverse and longitudinal) becomes 
of order one, meaning that, from now on, these fields can be also interpreted as incoherent superpositions 
of particles. The spectrum of these particles (mostly gluons) is obtained from the Fourier modes of the 
colour fields at time t > 1/Q S . These particles can interact with each other, as their density is still quite 
high (albeit decreasing with time, due to expansion). As we shall see, these interactions are expected to 
lead to a phase of local thermal equilibrium — the quark-gluon plasma (QGP). 

The intermediate, non-equilibrium, form of matter, which interpolates between the CGC in the 
initial wavefunctions and the QGP at later stages is known as the glasma (a name coined as a combination 
of 'glass' and 'plasma') (70). The main, qualitative, feature of the glasma is the presence of longitudinal 
colour flux tubes with transverse area ~ 1/Q^, which are boost invariant : the colour fields depend upon 
the proper time r but not also upon the space-time rapidity rj s . The variables r and rj s , which can be 
used instead of t and z = x 3 in the forward LC, are defined as 

1 , t + z , z 



Vt 2 - z 2 , T) s = - In => tanh rj s = - . (46) 

2 t — z t 
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Fig. 25: Left: the longitudinal colour flux tubes which develop in between the remnants of the two nuclei in the 
early stages (r < 1/Q S ) of a A+A collision. Right: various choices of coordinates in the future light cone. 



Under a boost along the z axis, r is invariant while rj s is shifted by a constant. Lines of constant r and 



of constant r] s are shown in Fig. [25] (right). The fact that the glasma fields depend upon r but not upon 
rjs is a consequence of the symmetries of the collision, as encoded in the classical field equations ( [39] ). 
This is also consistent with the hypothesis of uniform longitudinal expansion, as originally formulated 
by Bjorken. Specifically, Bjorken has assumed that (i) after being produced at t ~ z ~ 0, the parti- 
cles undergo free longitudinal streaming, meaning that they keep a constant velocity along the z axis; 
accordingly, the particles that can be found at some later time t at point z are those with a longitudinal 
velocity v z = z/t; (ii) the distribution of the produced particles is uniform in v z . Together, (i+ii) 
imply that the distribution at time t is independent of z/t, hence of r] s . Note that this argument iden- 
tifies the momentum rapidity y (cf. Eq. ( fTT] )) of the produced particles with their space-time rapidity 
T] s : tanhy = v z = z/t = tanh?7 s . Hence, the boost invariance of the glasma fields implies that the 
distribution of the particles produced by the decay of these fields is independent of y. This is a generic 
feature of the particle production at the classical level, that is, on an event-by-event basis: the associated 
spectra are boost invariant. But the physical spectra, as obtained after averaging the classical results with 
the CGC weight functions of the incoming nuclei, cf. Eq. ( [38] ), are rapidity-dependent, because of the 
respective dependencies of the weight functions, as introduced by the quantum evolution with Y. 

An interesting consequence of the above considerations, which might be related to a remarkable 



phenomenon seen in the RHIC [ |7T| - [74| and the LHC data |75] - [77| and known as the ridge, refers to the 
rapidity dependence of the two particle correlation. The latter is defined as 



ai Pa±i Vbi Pal. ) 



dN 2 



dVad 2 p a ± dr] b d 2 p b ± 



dN 



driad 2 p a ± I \ d7] b d 2 p b ± 



dN 



(47) 



and what is generally plotted is the ratio 1Z between this correlation and is disconnected part (below, iV a 
denotes the number of particles of type a in a given bin in pseudo-rapidity and azimuthal angle), 



K = 



(N a N b )-{N a ) (N b ) 

(Na) (N b ) 



(48) 



as a function of the rapidity and the azimuthal separations between the two particles, A77 = r\ a — r\ b and 
A<p — 4> a — 4> b . Remarkably, the data for A+A collisions at both RHIC and the LHC show the existence 
of correlations which extend over a large rapidity interval A77 ~ 4 4- 8, but restricted to small azimuthal 
separations A</> ~ (see Fig. [26] left). This means that particles which propagate along very different 
directions with respect to the collision axis preserve nevertheless a common direction of motion in the 
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freeze out 
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Fig. 26: Left: the 'ridge' in the di-hadron correlations as measured byRHIC (STAR). Right: the causal structure 
of the di-hadron rapidity correlations; the 2 hadrons A and B propagating with different rapidities were causally 
connected in the space-time region where the respective backward light-cones with vertices on the freeze-out 
surface were overlapping with each other. 

transverse plane. By causality, such a correlation must have been produced at early times, when these 
particles — which rapidly separate from each other — were still causally connected (see the right panel 
of Fig. [26]). A simple estimate gives 

|Ary[ 

7"max — ^freeze— out ^ 2 i (49) 

for the latest time at which these particles could have been correlated. For a freeze-out time 

^"freeze — out ^ 

10 fm/c, and rapidity separations A 77 > 4, one sees that these correlations must have been generated 
before 1 fm/c. 

Long-range rapidity correlations are natural in the glasma picture, where all the spectra (in partic- 
ular, the 2-particle ones) are independent of y — at least, at the classical level. After averaging with the 
CGC weight functions, as in Eq. ( [38] ), the 2-particle spectrum acquires a dependence upon the average 
rapidity of the two particles (r/ a + %)/2, but not upon their difference A77. This last argument remains 
correct so long as one can neglect quantum corrections due to soft gluon emissions within the rapidity 
interval A 77, which in turn requires A 77 < l/a s . Concerning the azimuthal collimation of the ridge, 
this is not a consequence of the glasma — the particles produced via the decay of the classical fields are 
emitted isotropically in the transverse plane in the rest frame of the medium — , but can be generated 
via radial flow : the local fluid element has some transverse velocity v±, which introduces a bias in the 
azimuthal distribution of the particles produced by the decay of a same flux tube (the final correlation is 
peaked around the direction of v±). Flow phenomena will be discussed in more detail in the next section. 

For sufficiently high energy, the long-range rapidity correlations invoked for A+A collisions 
should be also present in the p+p collision. In that case, one expects no flow, as there are fewer pro- 
duced particles and the freeze-out time is shorter. Yet, a small ridge has been measured in p+p collisions 
by the CMS collaboration at the LHC |78| , but only in the high-multiplicity events (which are believed 
to be more central) and within a limited range in p± (from 1 to 3 GeV), which is in the ballpark of the 
proton saturation momentum at the LHC. In that case, the azimuthal collimation could be explained by 
an intrinsic angular correlation between the emission of two particles from the glasma field (791. 

Let us also note another possible consequence of the glasma flux tubes: the presence of lon- 
gitudinal electric and magnetic fields implies the existence of topologically non-trivial configurations, 
characterized by a large density of Chern-Simons topological charge. Such configurations are interesting 
in that they break the charge-parity (CP) symmetry: via the chiral anomaly, they generate a difference 
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Fig. 27: The double-peak structure characteristic of elliptic flow: two wide peaks at A</> = and Acj) = tt which 
extend over a wide interval A 77 in rapidity. The strong correlation peak visible at (A0, A 77) ~ (0,0) is associated 
with the fragmentation of the trigger jet. 



between the number of quarks with right-handed and respectively left-handed helicity. In the context 
of HIC, they may generate a new phenomenon, known as the chiral magnetic effect 180] : quarks with 
opposite helicities can be separated by the ultra strong magnetic fields (B ~ 10 18 Gauss) created in the 
peripheral ultrarelativistic A+A collisions, thus leading to a charge asymmetry between the two sides of 
the reaction plane. Since the direction of the magnetic field varies from one collision to another, this 
effect leads to fluctuations in the distribution of the electric charge of the final hadrons. These theoretical 
expectations appear to be supported by measurements at RHIC (STAR) |8T|. 



4 The Quark Gluon Plasma 

The main topic of this chapter is the quark-gluon plasma (QGP) — the partonic form of QCD matter in 
thermal equilibrium which exists for sufficiently large temperatures, as demonstrated by numerical cal- 
culations in lattice QCD. This form of matter is expected to be created during the intermediate stages of a 
ultrarelativistic HIC, albeit only for a short lapse of time and in a state of only local thermal equilibrium. 
We shall first review the main experimental evidence in favour of QGP in a HIC, namely the observation 
of flow in the particle production and its successful description in terms of hydrodynamics. Then we 
shall discuss the QGP thermodynamics from the viewpoints of lattice theory and perturbative QCD. We 
shall also mention the difficulty of perturbation theory to describe the dynamics out-of-equilibrium (in 
particular, the transport coefficients and the process of thermalization). Then we shall consider the phe- 
nomenon of jet quenching (the energy loss by an energetic parton via interactions in the plasma), which 
is an important tool for exploring the deconfined matter produced in HIC's. At several places in what 
follows, we shall encounter situations where perturbation theory appears to be insufficient and which 
may signal a regime of strong coupling. To address such situations from the opposite limit — that of a 
coupling which is arbitrarily strong — , one can rely on techniques borrowed from string theory, via the 
AdS/CFT correspondence. This will be briefly discussed (in relation with the physics of HIC's) in the 
last section of these lectures. 



4.1 Correlations and flow in HIC 



In our previous discussion of the 'ridge', in Section we focused on the long-range rapidity corre- 
lations which are rather strongly peaked in Acf) = <j) a — <j) b near Acf) = 0, leading to the ressemblance 
with a mountain ridge. But as a matter of fact, the di-hadron correlations measured in A+A collisions at 
RHIC and the LHC show an even more pronounced double-peak structure, visible in Fig. [27j with large 
but wider peaks at Acf) = and Acf) = tt. (In the analysis leading to Fig. |26l this structure has been 
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Fig. 28: Left: a frontal view of a HIC, illustrating the impact parameter bj_, the almond-shape interaction region, 
and the azimuthal angle <p. Central and right: the distribution of nucleons within the interaction region in a given 
event can include elliptic (y<i), triangular (vs), or even higher harmonic modes. The reference angles fa, ^3 show 



the tilt of the interaction region with respect to the geometrical 'reaction plane' (cf. Fig. 23). 



subtracted away to render the ridge effect visible.) This is known as elliptic (or 'transverse') flow |82). 
As also visible in Fig. [27} this double peak structure can be well parameterized as 

dN P airs \ oc ^cos(2A0), (50) 



dA0 

with v 2 the 'magnitude of the elliptic flow'. 

The explanation of this phenomenon turns out to be quite simple: it reflects the anisotropy of the 



interaction region — the almond-shape region where the two nuclei overlap with each other; see Fig. 23 



for non-central collisions. This anisotropy entails a pressure gradient in the initial conditions: the 



pressure is larger along the minor axis of the ellipse (the x axis in the left panel of Fig. [28]) rather than 
along the major one; accordingly, more particles will be emitted in the direction of the largest gradient. 
This ultimately generates an anisotropy in the azimuthal distribution of the produced particles, which for 
symmetry reasons is of the form shown in Eq. ( [50] ): 

dN 

— - oc 1 + 2v 2 cos 2(f). (51) 
dcp 

This argument looks simple, as anticipated, but there is something deep about it: the role played by col- 
lective phenomena like pressure gradients or flow. Such phenomena are natural for many-body systems 
which relatively strong interactions — sufficiently strong to be able to transmit the asymmetry of the 
initial geometry into properties of the final state. This is an important point to which we shall return. 

The qualitative arguments above suggests that the coefficient v 2 characterizing the strength of the 
anisotropy should increase with centrality. This trend is indeed seen in the data (at least, for not too 
peripheral collisions, for which the interaction region becomes tiny and dilute). Also, particles which 
experience a stronger pressure gradient are expected to have a larger transverse momentum, as they 
inherit the velocity of the fluid; so, v 2 should rise with p±. This expectation, too, is confirmed by the 
data, at least for not too large p± < 5 GeV: very hard particles cannot be driven by the medium, so for 
them v 2 is naturally small. The measurements of v 2 (say, via 2-hadron correlations) yield very similar 
results at RHIC and the LHC (see the left panel of Fig. [30]), showing that v 2 is roughly independent of 
the COM energy. Note also that the typical values of v 2 for semi-hard momenta are relatively large, 
v 2 ~ 0.2, meaning that the collective phenomena alluded to above are indeed quite strong. 

But the elliptic flow and the ridge are not the only collective phenomena hidden in the di-hadron 



correlations illustrated in Fig. 27 By looking at the most central collisions where v 2 is relatively small, 
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Fig. 29: The 'ridge' at A<j) ~ and the 'double-hump' at \Acj) — ir\ ~ 1.1 as visible in the LHC data for di- 
hadron correlations in the 5% most central collisions (left panel), together with the harmonic decomposition of 
these correlations, cf. Eq. |52] ) (right panel) J#4|/. 



one sees not only the narrow 'ridge' at Acj) ~ 0, but also a 'double-hump' on the away side, at | Ac/)— tt\ ~ 



1.1, which extends too over a large interval At/ (see Fig. [29] left). The harmonic decomposition of this 
signal reveals higher Fourier modes with significant strengths, as illustrated in the right panel of Fig. [29 
This leads to the following generalization of Eq. ( |5T| ) : 

d ^-\ cxl + 2f>2)cos(™A0 (52) 

' 71=1 

where the various coefficients i; n up to vq have been extracted from the LHC data and they are compared 



to V2 in the right panel of Fig. [30] (All these coefficients are roughly independent of 77, meaning that they 
describe correlations over a wide interval A77.) 

What is the physics of such higher harmonics ? It is generally believed that they are the con- 
sequence of fluctuations in the distribution of nucleons within the interaction region, as illustrated in 
Fig. [28] (right panel) (83j. Namely, even though the overlapping region between the two nuclei has an 
elliptic shape, the nuclear matter inside it is neither homogeneous, nor strictly ellipsoidal, because of 
fluctuations in the particle distribution. The azimuthal distribution in a given event can be decomposed 
into harmonics, with coefficients v n and reference angles i\) n : 

dN 00 

— — oc 1 + y^2v n cosn(()) - ip n ) . (53) 
dcp ^ 

71=1 

The reference angles ip n are generally different from the conventional reaction plane (iJjrp = 0) and are 



difficult to measure, but they drop out in the 2-particle correlations, as manifest in Eq. ( [52] ). 

The fact that the initial geometry of the interaction region can have complicated fluctuations is not 
necessarily a surprise, given the granularity of the nucleons. What is remarkable though is the ability 
of the system to transmit these fluctuations into the distribution of the produced particles, via transverse 
flow. The effective theory for flow is hydrodynamics and will be succinctly discussed in the next section. 

4.2 Hydrodynamics and kinetic theory 

Hydrodynamics is the theory which describes the flow of a fluid independently of its detailed micro- 
scopic structure. More precisely, the equations of hydrodynamics have an universal form (at least, for 
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Fig. 30: Left: the RHIC and LHC results for as a function of p±, for collisions within the centrality bin 
40% - 50%; from Ref. [85]. Right: the LHC results (ATLAS) for v n (p ± ), for n = 2,3,4,5 and 6, and the 
centrality bin 20% - 30% JM. 



a given underlying fundamental theory, like QCD), but they involve a few 'parameters' which depend 
upon the nature of the fluid and can in principle be computed via microscopic calculations. The scope 
of hydrodynamics can be most easily explained with reference to thermodynamics. The latter describes 
a many-body system in global thermal equilibrium, in which the intensive quantities, like tempera- 
ture, pressure and the chemical potentials associated with the various conserved charges (electric charge, 
baryonic charge, etc) are time-independent and uniform throughout the volume V of the system. Hy- 
drodynamics can be viewed as a generalization of this picture towards a state of local equilibrium: the 
intensive quantities alluded to above can vary in space and time, but they do that so slowly that one can 
still assume thermal equilibrium to hold locally, in the vicinity of any point. Gradients of pressure and 
thermodynamics naturally lead to flow, with a local fluid velocity v which is itself slowly varying in 
space and time. 

The equations of hydrodynamics are simply the ensemble of the relevant conservation laws — for 
the energy, momentum and the other conserved charges: 

<9 M T^ = 0, M J£ = O, ••• (54) 

where T^ v is the energy-momentum tensor, is the density of the baryonic current (the volume integral 
of Jg is the difference between the number of baryons and the number of antibaryons), and the dots stand 
for other conserved charges. These densities depend upon the intensive (local) quantities describing the 
state of the fluid: the energy density e = E/V, the pressure P, the 4-velocity = 7(1, v) (with 7 = 
1/y/l — v 2 ), and a set of 'friction coefficients' known as viscosities, which characterize the dissipative 
properties of the medium. 

The relations between the densities of the conserved charges (T^, ...) and the intensive 
quantities are obtained via a gradient expansion with respect to the slow space-time variations of the 
latter. More precisely, this amounts to an expansion in powers of£/R, where R is a characteristic size of 
the system (in a HIC, R is the transverse size of the interaction region) and £ is the mean free path of the 
particles composing the fluid (the typical distance between two successive collisions). This quantity will 
play an important role in what follows, so let us open here a parenthesis and discuss it in more detail. 

At least for sufficiently weak coupling, the mean free path £ can be estimated using kinetic theory. 
This is an effective theory too, but it applies at shorter, microscopic, scales: in that context, the mean free 
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path is typically the largest scale in the problem (it is much larger than the Compton wavelength A ~ 1/k 
of a particle or the typical duration of a scattering processes). Kinetic theory allows one to follow the 
evolution of the particle distributions in phase-space — i.e. in space-time and in momentum space. 
To that aim, this theory involves more information about the microscopic dynamics, like cross-sections 
for the particles interactions described via the 'collision term' in the Boltzmann equation — the central 
equation of kinetic theory. But even without solving that equation, one can deduce an estimate for £ via 
simple considerations: the collision rate (the inverse of the typical time r co ii between two successive 
collisions) scales like r^Jj ~ nv Te \a, where n is the particle density, v Te \ is their average relative velocity, 
and a is the cross-section for their mutual interactions. The mean free path is then obtained as 

v 1 

^ ~ ^coii ~ ~ — , (55) 

v re \n a na 

where v is the average velocity of the particles, so v/v re \ is a number of order one. Since a is naturally 



proportional to some power of the coupling constant, Eq. ( |55| ) shows that the mean free path becomes 
smaller — meaning that the hydrodynamical description works better — when the coupling is strong. 

To be more specific, consider a system that will play an important role in what follows: a weakly 
coupled quark-gluon plasma with (local) temperature T. This is a nearly ideal gas of ultrarelativistic 
particles, so the particle densities scale like n ~ T 3 separately for quarks and gluons. To leading order 
in a s , scattering is controlled by the 2 — » 2 elastic collisions shown in Fig. [31} where the external lines 
represent thermal particles with typical energies and momenta of order T. These processes yield aoc^. 
However, for the processes involving the exchange of a gluon in the t channel, there is a logarithmic 
enhancement associated with the singularity of the Coulomb scattering at small angles: the Rutherford 
formula reads dcr/dQ oc o? s j(T 2 sin 4 9), with 9 the scattering angle, and it is strongly divergent when 
9^0. The cross-section a which is relevant for computing the mean free path ( [55] ) is not the total 
cross-section cr t ot = I dfHdcr/dfi), but rather the transport cross-section : 



= / <*(» — ») m K I M sln "< 1 " cos *> T&e ~ f* / T ~ ln i • <56) 



which more properly characterizes the efficiency of the interactions in redistributing energy and mo- 
mentum. The factor 1 — cos 9, which vanishes as 9 2 /2 at small angles, accounts for the fact that the 
small-angle scattering is inefficient in that sense, as intuitive from the fact that one cannot equilibrate 
an anisotropic energy-momentum distribution via collinear scattering. Due to this factor, the integral in 
Eq. ( [56] ) is only logarithmically divergent as 9 — » (unlike a to t, which would be quadratically diver- 
gent). In reality, this divergence is screened by plasma effects which occur at the momentum scale gT 



(see the discussion in Section 4.3). This implies that the minimal collision angles are 9 ~ gT jT ~ g, 
corresponding to transferred momenta of order gT, Hence, to leading logarithmic accuracy, the integral 
can be estimated as shown in the r.h.s. of Eq. ( [56] ). In turn, this implies the following estimate for the 
mean free path in a weakly-coupled QGP (cf. Eq. ([55])) 

T a z s m(l/a s ) 

which as long as g <C 1 is indeed much larger than both the Compton wavelength A ~ 1/T of the 
thermal particles and the typical duration ~ 1/gT of a scattering process. 

We now close the parenthesis dedicated to the mean free path and return to the discussion of 
hydrodynamics. As already mentioned, this is a legitimate effective theory for flow when £ <C R. The 
constitutive relation allows one to relate the energy-momentum tensor to the velocity field and its 
gradients, via an expansion in powers of£/R. The powers of £ are associated with dissipative phenomena, 
while those of 1 /R with the gradients in the fluid. To zeroth order in this gradient expansion one obtains 
the ideal hydrodynamics. This is 'ideal' in the sense that there is no dissipation. The corresponding 
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Fig. 31: 2—^2 partonic processes dominating the cross-section for (elastic) scattering at weak coupling. 



structure of T^ v follows entirely from the assumption of local thermal equilibrium. Namely, in the 
local rest frame of a fluid element (u^ F = (1, 0, 0, 0)), the energy-momentum tensor has the diagonal 
structure familiar from the thermodynamics: T£ F = diag(£, P, P, P). Boosting to the laboratory frame, 
where the fluid 4- velocity is u^, this yields 

T» v = (e + P)u^u v - Pg^ , (58) 



where g^ v = diag(l, —1, —1, —1) is the Minkowski metric tensor. The r.h.s. of Eq. (58) involves 5 



independent quantities: e, P, and the 3 components v % of the (local) velocity. Their space-time evolution 



is determined by the respective conservation law in Eq. ( [54| ), which yields 4 equations, together with the 
assumed equation of state, which specifies the functional relation e(P) between the energy density and 
the pressure (e.g., e = 3P for an ideal gas of massless particles). 

Since ideal hydrodynamics ignores dissipation, one might think that it corresponds to a situation 
where the coupling is weak, but that would be wrong: it rather corresponds to strong coupling. This may 
seem counterintuitive but it can be understood as follows: the dissipative phenomena are proportional to 
the ability of the system to transfer momentum in a direction perpendicular to the fluid velocity (since 
such a transfer results in slowing down the flow). Within kinetic theory at least, this transfer is realized 
by particles moving throughout the fluid in between successive collisions. Hence, the rate for transfer 
is proportional to the mean free path ( [55] ) and thus to the inverse of the coupling. As an example, 
consider the shear viscosity: this characterizes the friction force between two neighboring layers of fluid 
which propagate, say, along the x axis, but at slightly different velocities (so there is a non-zero gradient 
du x /dy). There is friction because some longitudinal momentum p x gets transferred from the faster 
layer to the slower one, at a rate proportional to the velocity gradient: 



1 dp x du x 

i^ = -^' (59) 



where A is the contact area between the two layers and 77 is the shear viscosity. Within kinetic theory, 
77 ~ Ipv ~ (p/n)(v/cr) where p is the mass density in the fluid and the second estimate follows after 
using Eq. ( [55] ). For a non-relativistic fluid p/n = m, so 77 ~ mv/cr, whereas for the weakly-coupled 
QGP, p = e ~ 3nT and v = 1 and therefore 

rp ^3 

V ~ ~ ~ 2 Wi / \ ' (60) 

a a^m{l/a s ) 

In both cases, the fluid density has canceled in the ratio p/n, so the viscosity is independent of the 
density, or, equivalently, of the pressure. (Recall that P = nfcgT, with ks the Boltzmann constant.) This 
remarkable conclusion has been first derived by Maxwell in 1860 via kinetic theory, and then confirmed 
by him experimentally. 



The l.h.s. of Eq. (59) represents a contribution to the component T xy of the energy-momentum 
(or 'shear') tensor : the flux of the x component of the momentum vector across a surface with constant 
y. So, Eq. (59) displays a dissipative correction to ; as expected, this is of linear order in the gradient 



expansion and it scales like £/R (since 77 ~ i and d y ~ 1/P). While first-order gradient corrections 
(leading to the Navier-Stokes equation) are sufficient to describe dissipation for a non-relativistic fluid, 
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this is not true anymore for a relativistic one: to be consistent with causality and Lorentz invariance, one 
must use a second-order formalism, which also includes quadratic terms in the gradient expansion. 

To summarize, for the problem of hydrodynamics to be well defined, one needs to specify (i) the 
equation of state e(P) (this is generally taken from lattice QCD calculations; see below), (ii) the time 
to at which hydrodynamical evolution can be turned on (meaning that local thermal equilibrium has been 
reached), (iii) the initial conditions at to for the energy density e{x) and the velocity v(x) fields, and 
(iv) the various viscosities like 77 which characterize the dissipative properties of the medium. Note that, 
in this context, the 'initial time' to is not the same as the time r s ~ l/Q s at which the CGC formalism 
provides the 'initial conditions' (in the sense of the discussion in Section [3J] ), but it is the a priori larger 
equilibration time r eq . Within the hydro simulations, this is a free parameter, like the viscosities or the 
parameters which enter the equation of state. These parameters are fixed a posteriori, by matching the 
results of the hydro evolution at the time of freeze-out onto some of the experimental results for particle 
production, like the centrality dependence of the particle multiplicities and of their average transverse 
momentum. 

For quite some time, roughly until 2007, it seemed that the RHIC data can be well accounted 
for (within the error bars) by ideal hydrodynamics [86]. This led to the conjecture that the deconfined 
matter produced at the intermediate stages of a HIC might be strongly interacting ('strongly coupled 
quark-gluon plasma' or sQGP). In order to test this conjecture, and also to describe the more accurate, 
recent data, it became necessary to include dissipative effects, within the second-order formalism. Full 



calculations in that sense, including comparison with RHIC data, became available only recently |87] - [89 | 
(and refs. therein). They are all consistent with a non-zero, albeit small, relative value of the viscosity, 
as measured by the ratio 77/5. Here, s is the entropy density, and the ratio 77/5 is dimensionless in natural 
units (in general, it has the dimension of h). This ratio is a natural measure of the deviations from ideal 
hydro, as we explain now. The entropy density s is proportional to the particle density; e.g., s = 4n for 
an ideal gas of massless particles. Thus, 77/5 ~ £v(p/n) ~ tj\, where A is the Compton wavelength 
of a particle in the fluid: A = l/(mv) in the non-relativistic case and A ~ 1/T for a weakly coupled 
QGR By the uncertainty principle, the ratio £/X cannot be smaller than h times a number of (9(1). So, 
the ratio 77/5 cannot become arbitrary small, even when increasing the coupling. In that sense, a physical 
fluid can never be ideal. 

An additional argument in that sense comes from the study of a strongly-coupled theory via the 
AdS/CFT correspondence. At least for the more symmetric, conformal, field theories to which it applies, 



this formalism predicts a lower bound on the ratio 77/5, namely |90,91| (the subscript 'CFT' refers to a 



conformal field theory; see Section 4.5 for details) 



> £ , (61) 

CFT 47r 



with the lower bound being reached in the limit of an infinitely strong coupling (in a sense to be char- 



acterized in Section [43] ). One remarkable thing about the heavy-ion data at RHIC and the LHC is that 
they seem to require a value 77/5 almost as small as this absolute lower bound: 77/5 ~ 0.08 0.20 
depending upon the details of the 'initial conditions' at time to. (For instance, the analysis in Ref. [[87) 
favors a value 77/5 ~ 0.16 for initial conditions of the CGC type, as shown in the left panel of Fig. [32]) 
The other remarkable thing is that, in order to be successful, the hydro descriptions of the data must 
assume a very small equilibration time to < 1 GeV/c. These are both hallmarks of a system with strong 
interactions. Indeed, the particles thermalize by exchanging energy and momentum (and other quantum 
numbers) with each other, via their mutual collisions. So, we expect the thermalization time to be shorter 
for strongly interacting systems. This expectation is supported by kinetic theory, which yields a thermal- 
ization time T eq ~ ijv oc \Ta 2 s ln(l/a s )] -1 to leading-order at weak coupling. For realistic values of 
a s , this perturbative estimate is too large to be consistent with the data (even when corrected for inelastic 
processes like 2—^3, which turn out to be important (92)). 
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Fig. 32: Left panel: the V2 results of hydro calculations using the second order formalism with CGC initial 
conditions and various values ofrj/s [87]. The comparaison with the RHIC data favors rj/s ~ 0.16, which is 
about twice the lower limit ( [6T] ) predicted by AdS/CFT at infinitely strong coupling. Right label: the longitudinal 
expansion has the tendency to collimate nearby particles, thus opposing to the evolution towards isotropy. 



In the context of heavy ion collisions, the evolution towards (local) thermal equilibrium is fur- 
thermore hindered by the extreme anisotropy of the initial conditions and also by the anisotropy of the 
early-time expansion, which is predominantly longitudinal. Recall the glasma picture of the initial condi- 



tions (at times r s ~ 1/ Q s ), which is that of colour flux tubes extending along the collision x, cf . Fig. 25 



Flux tubes have an internal tension opposing to their longitudinal extension, like a string. Accordingly, 
the longitudinal component of the energy-momentum tensor associated with the glasma is negative ; 
one finds T^ gma = diag(s, £, e, —e) at such early times, which is extremely anisotropic, as anticipated. 
Subsequent interactions among particles are supposed to restore isotropy, but this is rendered difficult by 
the longitudinal expansion, as illustrated in the right panel of Fig. [32] Namely, even if the particle distri- 
bution turns out to be locally isotropic at a given position in space and time, the subsequent anisotropic 
expansion rapidly separates the particles from each other according to the directions of their velocities: 
only those particles remain close to each other which had nearly parallel velocities. In other terms, by 
itself, the longitudinal expansion would naturally build a particle distribution in which nearby particles 
move along quasi-parallel directions, thus opposing isotropy. To beat this tendency and ensure isotropy, 
one needs strong interactions which continuously randomize the directions of motion of the particles. 
This would be natural at strong coupling, as alluded to above. But there are also other scenarios which 
are currently explored, including weak-coupling ones. One promising mechanism in that sense refers to 
plasma (Weibel) instabilities : due to the anisotropy of the expanding parton distribution, the soft colour 
fields radiated by these partons can develop unstable modes, that is, modes whose amplitudes grow expo- 
nentially with time (at least, during a limited time interval). So far, it is not clear whether this mechanism 
can lead to rapid isotropisation in the presence of longitudinal expansion, but its studies are under way 
(see Refs. (93 -98| for recent work and related references). Recent developments include a calculation 
(similar to previous work in inflationary dynamics) of the spectrum of initial quantum fluctuations in the 
glasma (95), a parametric analysis of the interplay between plasma instabilities and Bjorken expansion in 
the weak-coupling limit (96}|97J, and an interesting scenario (still at weak coupling) in which the elastic 
scattering between the highly occupied glasma fields leads to the formation of a transient Bose-Einstein 
condensate (99). 
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Fig. 33: Schematic representation of the phase-diagram in QCD at finite temperature and non-zero quark density, 
as emerging from lattice calculations at zero (or small) quark density and from various theoretical considerations 
(like pQCD) in the other domains. 



4.3 QGP: Thermodynamics and collective excitations 

In this section, we shall deviate from the experimental situation in HIC's, where the partonic medium is 
rapidly expanding, and focus on a quark-gluon plasma at rest, in thermal and chemical equilibrium. The 
existence of such a deconfined phase in QCD at finite temperature has been unambiguously demonstrated 
via numerical calculations on a lattice, which have also given a lot of information about the thermody- 
namics of this system. Some of this information has been corroborated via analytic calculations at weak 
coupling, which turned out to be very non-trivial. The analytic methods become essential when one 
is interested in real-time phenomena, like the response of the system to time-dependent external per- 
turbations, as characterized by transport coefficients. Indeed, real-time phenomena cannot be (easily) 
studied via lattice calculation^} which are a priori formulated in a space-time with Euclidean signature 
('imaginary time'). It should be also stressed that, even though weak-coupling techniques appear to be 
quite successful in reproducing the lattice results for the QGP thermodynamics, the hypothesis that the 
coupling be strong is not yet totally excluded (within the temperature range relevant for the phenomenol- 
ogy at RHIC and the LHC): indeed, weak-coupling calculations seem unable to explain the small r\j s 
ratio supported by the data (cf. Section |4~2| ). In what follows, all that will be explained in some detail. 

Fig. [33] shows a cartoon of the phase-diagram expected in QCD when varying the temperature 
T and the net quark density (or the quark chemical potentials by which one means the difference 



between the density of quarks and that of the antiquarks. (For simplicity, Fig. [33] treats the three light 
quark flavors — the only ones to be relevant for the phase diagram — on the same footing.) This 
diagram has been actually demonstrated only in special corners, like the deconfinement phase transition 
with increasing T at zero (or small) fermionic density, that has been established on the lattice, and the 
islands denoted as 'nuclei' or 'neutron stars', which are rather well understood within nuclear theory. 
The 'colour superconductivity' phase at high quark density, which is predicted by pQCD (at least for 
p » AqcdX will n °t be discussed here, as it is not expected to play any role in the ultrarelativistic 



6 There is some recent progress in computing transport coefficients on the lattice (see e.g. the review paper 1 100] and Refs. 
therein), but although promising, this method is still inaccurate and very fastidious. 
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Fig. 34: The quark-antiquark potential between two heavy quarks, as computed in QCD on the lattice. Left panel: 
T = 0. Right panel: various temperatures, which are all larger than the critical temperature for deconfinement T c 
(for comparison, the T = potential is also shown, as the continuous line which keeps rising) f 1 04] /. 



HIC's. (See the review papers ]101||102[ and Refs. therein for detailed discussions of this phase.) 

As also illustrated in Fig. [33] the deconfinement phase transition has been first explored during the 
expansion of the Early Universe: the high temperature 'soup' of matter created right after the Big Bang 
was originally in the deconfined, QGP, phase; due to its rapid expansion, this matter has cooled down 
and thus crossed into the confined, hadronic, phase, at a very short time ~ 10 -5 seconds after the Big 
Bang. In the context of HIC's, this transition is being probed the other way around: to start with, the 
partons are confined within the nucleons composing the two nuclei; the collision liberates these partons 
and, if their energy density is high enough, they can thermalize at a temperature superior to the critical 
temperature for deconfinement. If so, they form a transient QGP phase which cools down via expansion 
and eventually 'evaporates' into hadrons. In both scenarios, the net quark density is small and plays no 
role for the transition. In the Early Universe, the excess in the number of quarks over antiquarks was 
negligible (if any !) [ 103 ]. In HIC's, there is of course a net baryon number, due to the 2 A ~ 400 
nucleons within the incoming nuclei; however, this excess is small compared to the number of hadrons 
(a few thousand) produced in the final state. This implies that most of the partons which exist in the 
intermediate stages of the collision are actually gluons or 'sea' quark-antiquark pairs. 

The fundamental property of the QGP is, of course, deconfinement : quarks and gluons can move 
(more or less freely, depending upon their mutual interactions) throughout the whole volume of the 
plasma, without being confined within hadrons with radia ~ 1/Aqcd- How is this possible ? A quark 
and an antiquark in isolation (i.e. at zero temperature) attract each other via a force which becomes 
roughly constant — corresponding to a linear potential; see Fig. [34] left — at distances r > 1/Aqcd- 
Due to this force, the qq pair is tightly bound ('confined') into a meson. This meson can be broken, 
say, via a hard scattering, but only at the expense of producing additional gluons and qq pairs which 
'glue' to the original quark and antiquark, in such a way to form colour singlet states (new hadrons). 
This is the situation in the 'usual' hadronic processes, including p+p collisions at the LHC, where the 
parton density right after the collision is not very high — so, these partons can evolve and eventually 
hadronize independently of each other. But in HIC's, the density of the liberated partons is such that the 
typical interparticle separation is much shorter than 1/Aqcd- At such short distances, the attraction force 
between these partons is smaller than their kinetic energy, so the partons can move around each other and 
arrange themselves in such a way to minimize their mutual repulsion. The net result is that the colour 
charge gets screened over relatively short distances r <C 1/Aqcd> thus preventing the development of 
confining forces at larger distances. 
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Fig. 35: Left panel: illustration of the Debye screening in a QED plasma; the screening length is the inverse of 
the Debye 'mass': td = 1/rriD- A similar mechanism is active in the QGR Right panel: the Debye screening in 
QCD is the result of one-loop corrections due to the 'hard' (k ~ T) thermal particles (quark and gluons) which are 
resummed in the propagator of the 'soft' (k ~ gT) electric gluon exchanged between the quark and the antiquark. 



This colour screening in the QCD plasma is very similar to electric (Debye) screening in ordinary, 
electromagnetic, plasmas, or in electrolytes. Ions with positive electric charge attract ions or electrons 
with negative charge, in such a way to form clouds of particles which look electrically neutral when 
seen from far away: the net charge decreases exponentially with the distance from the central charge 



(see Fig. [35] left). In the context of QCD, the 'positive and negative electric charges' are replaced by the 
Af? — 1 = 8 'colour' charges carried by quarks and gluons, but the exponential screening of the chromo- 
electric charges works in a similar way. As a consequence of that, the quark-antiquark potential flattens 
out (meaning that there is no attraction force) at large distances r > with rr> the Debye radius. This 
flattening is clearly seen in the lattice calculations at finite temperature (see e.g. the discussion in | |104[ ), 



as illustrated in the right panel of Fig. [34] This also suggests that in a finite temperature plasma one 
cannot have quarkonia (bound states made with a heavy quark (Q) and a heavy antiquark (Q), like J/^, 
with size Tqq > rrj. This observation [ 105J led to the fertile idea of quarkonia melting in a quark-gluon 



plasma, a very active field of research for both theoretical (including lattice) and experimental studies of 
HIC's. (See [106] for a recent overview of the theory and more references.) 

The lattice calculations also allow one to study the deconfinement phase transition with increasing 
temperature. The respective results for the pressure and energy density are illustrated in Fig. [36] They 
show a sudden increase around a critical temperature T c ~ 160-1-180 MeV, interpreted as the result of the 
rise in the number of degrees of freedom (d.o.f.), due to the liberation of quarks and gluons. For T < T c , 
the only 'thermodynamically active', hadronic d.o.f. (those whose masses are not much higher than T) 
are the 3 pions: tt° and tt ± . For T > T c , this number jumps from 3 to 52: the gluons, which appear in 8 
colours and 2 transverse polarizations (8x2 = 16 d.o.f.), and the 3 light quarks and antiquarks, each of 
them having 2 spin states and 3 possible colours (3x3x2x2 = 36 d.o.f.). Lattice calculations become 
more tedious for light quark masses and the extrapolation to physical quark masses has become possible 
only recently [ |107[|108[ . This is important since both the actual value of the critical temperature and the 
nature of the phase transition are strongly influenced by the values of these masses. A phase transition 
is said to be of nth order if it involves a discontinuity in the derivative of order (n — 1) of the pressure. 
For instance, if the QCD phase transition was of first-order, then it would proceed via a mixed phase 
where hadronic bubbles coexist with regions of QGR But this is not what happens in QCD: recent lattice 



calculations [109| show that, for physical quark masses, the deconfinement phase transition is truly a 
cross-over, that is, a relatively smooth process during which the pressure and all its derivatives remain 
continuous across the transition. 

But albeit smooth, the phase transition represents a genuinely non-perturbative phenomenon, 
which cannot be described within perturbative QCD. This can be appreciated e.g. by inspection of 
the lattice results for the trace anomaly (e — 3P)/T 4 , which exhibit a sharp peak around T c , as visible 
in the left panel of Fig. [37] As mentioned in Section [4~2] e = 3P for an ideal gas of massless particles, 
meaning that the corresponding energy-momentum tensor = diag(^, P, P, P) is traceless: = 0. 
This property is in fact a general consequence of conformal symmetry : it holds for any theory which 
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Fig. 36: The lattice results (as obtained by the Budapest-Wuppertal collaboration JJlOj /j for the QCD pressure and 
energy density as a function of the temperature. 



involves no intrinsic mass parameter and hence is invariant under dilations. This is in particular the case 
for QCD with massless quarks at the classical level. However, at the quantum level, conformal sym- 
metry in QCD is broken by the radiative corrections responsible for the running of the coupling, which 
introduce the mass scale Aqcd- So, not surprisingly, the corresponding 'trace anomaly' (the deviation 
of from zero) is proportional to the /3-function, which measures the running of the coupling : 

dP _ x dg 



rpfj, 



3P 



(62) 



dg ' 1 <91n/x ren ' 

Here /i ren is the renormalization scale, as introduced by the subtraction of the ultraviolet divergences. In 
perturbation theory at finite temperature, it is convenient to choose fi ren = 2ttT as the central value and 
study the dependence of the results upon variations in ^ ren (typically by a factor of 2) around this central 
value. These variations measure the stability of the calculations against higher order corrections and thus 
are indicative of the theoretical uncertainties. They are shown as 'error bands' for the theoretical results 
in Figs.[37]and[39l 



The peak in the l.h.s. of Fig.[37]is a hallmark of the phase transition and is clearly non-perturbative. 
But for temperatures above T c , the 'trace anomaly' is rapidly decreasing (its relative strength becomes 
of order 10% for T > 3T C ), thus suggesting that a perturbative approach may become viable. This is 
furthermore indicated by the fact that, for T > 3T C , the pressure and the energy density reach about 80% 
of the respective values for an ideal gas of quarks and gluons, denoted as 'SB' (from 'Stefan-Boltzmann') 



in Fig. [36J A deviation of 20% may seem sufficiently small to be easily accommodated in perturbation 
theory, but this turns out not to be the case. There are two main reasons for that. First, unlike what 
happens at T = 0, where perturbation theory in QCD is an expansion in powers of a s = g 2 /4:7T, at finite 
temperature this is rather an expansion in powers of g, for reasons to be shortly explained. Second, the 
relevant values of the QCD running coupling are not that small: for T ~ 3T C and hence 2ttT ~ 2 GeV, 
one has a s ~ 0.25 and hence g = 1.5-1-2. For such large valus of g, there is no reason why an 
expansion in powers of g should converge, and indeed it does not: as visible in the right panel of Fig. [37} 
the successive corrections of 0(g 2 ), 0(g 3 ), 0(g 4 ), and 0(g 5 ) jump up and down, without any s ign o f 
convergence. (The weak-coupling expansion of the pressure in QCD is known to 0(g 6 ln(l/#)) [ 112| , 
which is the highest order that can be computed in perturbation theory: the corrections of 0(g 6 ) and 
higher are afflicted with severe infrared divergences due to magnetic gluons; see below.) 

The reason why, at finite temperature, perturbation theory is an expansion in powers of g rather 
than a s is because the quantum corrections associated with soft gluons — those with momenta k much 
smaller than T — are amplified by the Bose-Einstein thermal distribution function: 

1 T 

n B(k) = -= ~ — > 1 when E k = \k\ <C T. (63) 
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Fig. 37: Left panel: the trace anomaly in QCD, as numerically computed on the lattice {107\ 108]) , together with 
the respective predictions of a 'HTL-resummed' perturbation theory jlll]j (see the text for details), shown as a 
band (due to the uncertainty in the choice of the renormalization scale fjL ren ). Right panel: the predictions of (strict) 
perturbation theory for the pressure, as computed up to order g 2 , g 3 , g A , and g 5 , respectively. Whereas theO(g 2 )- 
result appears to match the lattice results (the dark grey band) rather closely, this agreement is spoilt after including 
higher order corrections; the ensuing series in powers of g shows no sign of convergence. 



This property is generic: it holds for any field theory which involves massless bosons (e.g., it holds 
for photons in a QED plasma). When k —> 0, the thermal factors nsik) — T/k lead to infrared 
divergences in the calculation of Feynman graphs, which are regulated by plasma effects, like Debye 
screening. These effects typically enter at the 'soft' scak^JgT. For instance, the Debye mass mp = 
which characterizes the exponential screening of the electric colour charge by the plasma constituents, 
is generated by the one-loop diagrams illustrated in the right panel of Fig. |35j which yield rnrj ~ gT. 
The resummation of these diagrams within the propagator of the exchanged gluon, as also illustrated in 

Fig. 



35 



: E k = Jk* + 



rn 



D' 



Hence, when k — >> 0, the 



renders the (electric) gluons effectively massive 
Bose-Einstein occupation number remains finite, but it is parametrically large: ng(fc) ~T/mrj ~ V#- 
This inverse power of g changes the perturbative order of the 2-loop correction to the pressure with one 
'hard' loop (k ~ T) and one 'soft' (k < gT) from a 2 s ~ g A to g A riB(k) ~ g 3 . This is the origin of the 
odd powers of g in the perturbative expansion. 



One should also mention here that Debye screening, as illustrated in the left panel of Fig. 35 is op 



erational for the electric gluons (i.e. for the Coulomb interactions), but not also for the magnetic ones — 
those having transverse polarizations. For non-relativistic plasmas, magnetic interactions are suppressed 
by powers of the velocities, but for (ultra)relativistic plasmas, like the QGP, they are as important as 
the electric ones. One expects magnetic interactions in the QGP to be screened at the 'ultrasoft' scale 
^mag ~ g 2 T, but the associated physics — in particular, the contribution of the 'ultrasoft' magnetic 
gluons to thermodynamics, which starts at 0(g 6 ) — cannot be computed in perturbation theory. Indeed, 
each additional 'ultrasoft' loop is accompanied by a factor ~ g 2 riB(k) ~ 5 f2 (7 1 /m ma g) ~ 1, meaning 
that diagrams with arbitrarily many such loops contribute at the same order in g. This explains why the 
corrections of 0(g 6 ) are non-perturbative, as alluded to above. 

In QCD, the Debye mass is only one example of a class of one-loop 'corrections' which are non- 



7 This scale gT is truly 'soft' only so long as g <C 1 ; as we shall shortly argue, results obtained under the assumption that 
the coupling is weak can be extrapolated towards g ~ 1 provided the plasma effects are properly taken into account. 

8 There is strictly speaking a difference between the Debye mass itid, which governs the infrared (k —> 0) limit of the static 
(ko = 0) propagator for the electric gluons, and the thermal mass m g , which enters the dispersion relation for the on-shell 
gluons; but these quantities are proportional with each other and are both of order gT; see e.g. Bl 13141 151 . 
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Fig. 38: A few examples of hard thermal loops (HTL's). All the external momenta are soft (ki ~ gT) and the loop 
integrations are dominated by relatively hard modes with p ~ T. 



perturbative at the 'soft' scale gT and should be viewed as a part of the leading-order theory at that 
scale, and not as corrections [ |113||116[ . These diagrams are known as hard thermal loops, since the 



typical momenta within the loop are of order T (the value preferred by the statistical, Bose-Einstein 
and Fermi-Dirac factors) and thus are hard compared to the soft (k ~ gT) momenta flowing along the 
external legs. There are HTL's with any number n of external gluons lines and with either zero, or two, 



quark external lines (see Fig.[38]for some examples). They are generally non-local, that is, they depend 
upon the external, soft, momenta. In particular, the HTL's for the 2-point functions (the quark and gluon 
self-energies) encode phenomena like Debye screening for the electric gluons, dynamical screening (or 
Landau damping) for the magnetic gluons, and the dispersion relations for on-shell quanta with momenta 
of order gT. Such quanta have wavelengths \ ~ 1/gT which are parametrically larger than the typical 
separation ~ 1/T between the typical plasma constituents — quarks and gluons with momenta of order 
T. Accordingly, the soft modes are truly collective excitations (or 'plasma waves'), with either quark or 
gluon quantum numbers. 

As already mentioned, the HTL's are of the same order as the respective tree-level amplitudes with 
'soft' external legs, so they cannot be expanded out in perturbation theory. Rather, one needs to perform 
a reorganization of the perturbation theory in which the HTL's are viewed as a part of the leading-order 
theory for the soft modes. Roughly speaking, this amounts to expanding around a gas of dressed quasi- 
particles whose zeroth-order properties (propagators and interaction vertices) are encoded in the HTL's. 
In practice, there are various ways to perform such reorganizations and it is quite reassuring that all the 



methods that have been proposed so far ]111||112|114||115||117||118J appear to be successful in describing 



the lattice data (although with considerably different amounts of efforts). 

One of these methods, known as 'HTL perturbation theory' (HTLpt) |118[ , consists in including 
the HTL's in the 'tree-level' effective theory, by adding and subtracting £htl (the sum of the HTL 
amplitudes) to/from the original Lagrangian: 

^QCD = A) + Ant = (A) + £rTl) + ( ~~ ^HTL + Ant) = £>Q + ^int ■ (^ 4 ) 

In this equation, Co is the free (g = 0) piece of the QCD Lagrangian, Ant is the respective interac- 
tion piece, C f = Co + £htl represents the new 'tree-level Lagrangian' which defines the Feynman 
rules (HTL-resummed propagators and vertices) of HTLpt, and, finally, C[ nt = Ant — £htl is the new 
'interaction Lagrangian'. The subtracted piece — £htl within C[ nt acts as a 'counterterm' to prevent 
double counting (the HTL's have been already included in the zeroth-order theory and they should not 
be regenerated via loop corrections in HTLpt) and also to correct for the fact that, within C f , the HTL's 
are used for all the modes, including the hard modes to which they do not really apply (this introduces 
spurious contributions at lower orders which are compensated by the 'counterterm' only in higher or- 
ders). In order for such compensations to efficiently work, one needs to go up to relatively high orders in 
HTLpt, which involve very tedious calculations (due to the non-local nature of the HTL's). It was only 
recently, after pushing such calculations up to three loop order that one has finally reached a good 



52 




£0.2 



0.0^ 



600 
T (MeV) 



Fig. 39: Comparison between the predictions of two versions of HTL-resummed perturbation theory for thermo- 
dynamics and the respective lattice results. Left panel: the predictions ofHTLpt fllff| / for the pressure at one-loop 
(LO), 2-loop (NLO) and, respectively 3-loop (NNLO) order fl ll]j vs. the lattice results (small cercles, triangles 
or squares) from two different collaborations. Right panel: the 2-loop result of the '2-particle irreducible resum- 
mation' of the entropy (11 7]j — solid and dotted lines correspond to two successive approximations for the thermal 
masses — vs. the lattice results shown as the grey band. In both cases, a good agreement with the lattice results 
is observed for temperatures T > 2.5T C . The theoretical 'error bands' follows from varying the renormalization 
scale in the range ttT < /i ren < AirT. 



agreement with the respective lattice results for T > 2.5T C . This agreement is visible in the left panel of 
Fig. [37] for the trace anomaly and in the left panel of Fig. [39] for the pressure. 

A more economical approach, in which a similarly good agreement with the lattice results (see the 
right panel of Fig. [39]) has been obtained via a simpler, 2-loop, calculation, is the '2-particle irreducible 
(2PI) resummation' of the entropy [ |114||117[ . In that approach, most of the difference between the 
2PI result for S and the corresponding result Ssb for the ideal gas comes from the thermal masses 
m qi m g ~ gT acquired by the hard (k ~ T) quarks and gluons via interactions in the plasma. Albeit 
formally small (m q ^ g ~ gT <C k ~ T), these masses cannot be expanded in perturbation theory, since 
such an expansion would generate powers of m 2 /k 2 leading to infrared divergences in the integral over 
k. The success of the 2PI description supports the physical picture of the QGP in terms of quasi-particles 
— quarks and gluons with typical momenta k ~ T, which are dressed by the medium (in particular, in 
the sense of acquiring thermal masses), but whose residual interactions are relatively small. Moreover, a 
substantial part of these residual interactions can be associated with collective excitations and screening 
effects at the 'soft' scale gT, as encoded in the HTL-resummed propagators. 

It is furthermore interesting to notice that the HTL resummation is based on the separation of 
scales gT «T which is a priori valid at weak coupling (g <C 1), yet this turns out to rather successfully 
describe the lattice data in a temperature range where g = 1.5-1-2. This confirms that the failure of 
ordinary perturbation theory (cf. Fig. [37| right panel) is not imputable to the fact that the coupling is 
relatively strong, but rather it is a consequence of expanding out the medium effects in powers of g in a 
kinematical domain where they are truly non-perturbative. 

Yet, the issue of the strength of the coupling remains open in so far as the study of dynamical 
phenomena is concerned. These phenomena refer to non-trivial evolutions in time, say off-equilibrium 
deviations in response to external perturbations. As long as the perturbations are small, their effects can 
be computed within the linear response theory, via the Kubo formula: the response of the plasma is linear 
in the strength of the perturbation, with a proportionality, or 'transport', coefficient which represents a 
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correlation function in thermal equilibrium. For instance, a constant electric field E acting on the quark 
constituents of the plasma (which carry electric charge) induces an electromagnetic current with density 
Oem) = &E l , where a is the electrical conductivity. The respective Kubo formula relates a to the long- 
wavelength (k l —¥ 0) and zero-frequency (u — » 0) limit of the current-current correlator in thermal 
equilibrium. Other transport coefficients include the (quark) flavor diffusion coefficients and the shear 
viscosity 77 introduced in Eq. ( [59] ). The use of Kubo formulae for perturbative calculations at weak 
coupling turns out to be quite tedious, because of the need to perform sophisticated resummations |119| . 



However, these formulae are very useful for non-perturbative calculations of the transport coefficients, 



either on the lattice 1 100], or via the AdS/CFT correspondence at strong coupling [120[ . 

For a weakly coupled QGP and to leading order in the coupling, the transport coefficient can be 
alternatively, and more efficiently, computed from the Boltzmann equation (linearized with respect to 
the off-equilibrium perturbation). This amounts to solving a linear integral equation which effectively 
resums an infinite number of diagrams of the ordinary perturbation theory in thermal equilibrium. These 
diagrams describe multiple scattering via soft gluon exchanges and can be generated by iterating the 
2 — 2 elastic processes shown in Fig. [31] arbitrarily many times. The ensuing transport coefficients are 
of the parametric form anticipated (on the example of the shear viscosity) in Eq. ( |60| ), but the use of 
the Boltzmann equation allows one to obtain more precise results, which are complete to leading order 



in a s [121[|122[ . Yet, these results are deceiving with respect to the heavy-ion phenomenology: as 



already mentioned in Section 4.2 the leading order estimate for 77/5 is too large to be consistent with the 
hydrodynamical description of the data. 

The last observation raises the question of the next-to-leading order corrections. Their calculation 
is extremely complicated and so far this has been accomplished for just one quantity: the diffusion 
coefficient D for a heavy quark with mass M » T. In the context of HIC, this quantity controls the 
collisional energy loss and the thermalization of heavy quarks like the charm or the bottom. Once again, 
the LO perturbative estimate for D [ |123[ appears to be too large to be consistent with the data. The NLO 



correction to D is of relative order g and has been computed in Ref. [124]. This appears to go in the right 
direction (it diminishes the value of D), but the effect is extremely large for realistic values of g — the 
NLO 'correction' is almost an order of magnitude larger than the respective LO result ! — , thus rising 
doubts about the reliability of the whole scheme. It looks like the perturbative series suffers from a lack- 
of-convergence problem similar to that noticed for the pressure. It might be that this problem too will be 
cured by all-orders resummations of the HTL's; but this issue is still open since such resummations have 
not yet been performed for dynamical quantities. Alternatively, there is the possibility that the transport 
phenomena, which involve long-range dynamics, be sensitive to rather large values of the QCD running 
coupling, which exclude weak-coupling techniques. If so, one could search for physical guidance in the 



corresponding results at strong coupling, as obtained via the AdS/CFT correspondence (see Section [43 
below). Finally, let us notice that the first lattice results for the transport coefficients have started to 
emerge, although the current errors bars are still quite large. These calculations are very difficult as they 
require to numerically perform an analytic continuation (from imaginary time to real time), which in 
turns requires very precise numerical data. In view of that, it is quite encouraging that the recent lattice 
results for the heavy quark diffusion coefficient [ |125[|126[ appear to be consistent with the heavy-ion 
phenomenology, within the (lattice and experimental) errors bars. 

4.4 Jet quenching 



In Section 3.6 we have mentioned two interesting phenomena occurring in 'dense-dilute' (p+A or d+A) 
collisions — the suppression of particle production and that of azimuthal di-hadron correlations at for- 
ward rapidities — , which in that context have been interpreted as consequences of gluon saturation in 
the wavefunction of the nuclear target: the larger the rapidity, the smaller the values of the longitudinal 
momentum fraction that are probed in the nucleus, and hence the stronger the saturation effects. On the 
other hand, the RHIC data for d+Au collisions at central rapidities (77 < 1) show no similar suppression 
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Fig. 40: The nuclear modification factor Raa, as measured in central Au+Au collisions at RHIC (PHENIX, left 
panel) and in Pb+Pb collisions at the LHC (CMS, right panel). For comparaison, one also shows the data for the 
production of direct photons, which show no nuclear effect as expected (left panel) and for hadron production in 
peripheral collisions, for which the nuclear effects are quite small (right panel). 



(see the 2 left-most plots in Fig. [18] and also the corresponding data in the right plot in Fig. |4T] below), 
which implies that nuclear saturation effects are not important in the central-rapidity kinematics at RHIC. 
(But this is likely to change at the LHC; see e.g. (14||69}[T27J.) The situation is however different for 
the 'dense-dense' A+A collisions: the respective data at RHIC and the LHC show a strong suppression 
of particle production and of the azimuthal di-hadron correlations already for central rapidities (and for 
relatively hard transverse momenta). These phenomena cannot be related to 'initial-state', saturation, 
effects in the incoming nuclei (at least not fully), since they do not show up in the mid-rapidity d+Au 
data at RHIC. Rather, they must correspond to interactions in the final state, that is, interactions with the 
dense partonic medium (the glasma and the quark-gluon plasma) which exists at intermediate stages. 
The change in the properties of a hard particle or of the associated jet induced by its interactions in the 
medium is generally referred to as jet quenching (see e.g. [128f[13l| for recent reviews). 

The suppression of particle production in A+A collisions is best characterized by the ratio Raa 



(the 'nuclear modification factor'), defined by analogy with Eq. (40), that is 



Raa(v,P±) = 



Ncoll 



dN h 




d 2 p±dr] 


AA 


dN h 




d 2 p_\_dr] 


pp 



(65) 



where the number N co u of binary collisions at a given impact parameter scales like A 4 / 3 (for relatively 
central collisions): indeed, there is a factor A 1 / 3 associated with the longitudinal width of each of the two 
nuclei and an additional factor of R 2 A oc A 2 / 3 coming from the integral over all the impact parameters. 
The experimental results for Raa (at mid-rapidities) in Au+Au collisions at RHIC (^/snn = 200 GeV) 
and in Pb+Pb collisions at the LHC (^/s^n = 2.76 TeV) are shown in the left and right panels of Fig. 40 
respectively. As anticipated, they show a substantial suppression of the hadron production as compared 
to p+p collisions, which persists up to p± ~ 20 GeV, at least. The interpretation of this suppression 
as a dense-medium effect is furthermore supported by the fact that (i) the direct photons (which do 
not interact with the hadronic matter) show indeed no suppression (cf. the left figure), and (ii) even 
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Fig. 41: Left: cartoon of a typical di-jet event as produced by a hard scattering. Right: azimuthal distribution of 
hadrons with p±_ > 2 GeV relative to a trigger hadron with p±_ > 4 GeV, as measured at RHIC (STAR), in p+p, 
d+Au and central Au+Au collisions [132]. 



for hadrons, the suppression is considerably smaller in the peripheral collisions (cf. the right figure), in 
agreement with the fact that the density and the size of the produced medium are much smaller in that 
setup. Note also that the suppression at intermediate values p± ~ 6 -=- 7 GeV is stronger at the LHC than 
at RHIC, indicating that the medium produced there is denser, as expected. 

Concerning the suppression of azimuthal correlations in the di-hadron production, this is clearly 



visible in the right plot in Fig. 41 which shows data taken at RHIC for hadrons with p± > 2 GeV : 
unlike for p+p and d+Au collisions, where one can see a peak at A<3> = tt, as expected for a pair of 
hadrons which are produced back-to-back, there is no such a peak in the central Au+Au collisions. This 
is interpreted as the consequence of the interactions suffered by the 'away' particle (the one that would 
have normally emerged at A$ = tt) while propagating through the medium. Via such interactions, the 
particle transverse momentum has been degraded and/or the particle has been deviated towards different 
directions, so it will not show up around A<3> = tt, nor in the original bin in p±. 
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Fig. 42: Left panel: a highly asymmetric di-jet event in Pb+Pb collisions at the LHC as measured by CMS. Right 
panel: a cartoon of an asymmetric di-jet event in A+A collisions. The hard scattering producing the jets occurs 
near the edge of the fireball. One of the jets (the 'trigger jet') leaves the medium soon after its formation and thus 
escapes unscattered, while the other one (the 'away jet') crosses the medium and is strongly modified by the latter. 
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Fig. 43: The distribution of the CMS data for p+p collisions (y/s = 7 GeV) and Pb+Pb collisions (y/s^N — 
2.76 TeV) as a function of Aj (left panel) and respectively of the azimuthal separation A<£i 2 w.r.t the leading-jet 
axis (right panel) in different bins of centrality: (a) refers to p+p, (f) to the 10% most central Pb+Pb, (e) to the 
centrality bin 10-20%, etc. The histograms show the respective predictions ofPYTHIA, which agree well with the 
data for p+p and for peripheral Pb+Pb, but not also for central Pb+Pb. From Ref. [134]. 



Besides confirming and sharpening the discoveries at RHIC, the first heavy ion data at the LHC 
revealed a new phenomenon, whose observation was possible because of the unprecedented ability of 
the detectors there (notably the calorimeters at ATLAS and CMS, with a wide coverage in rapidity) to 
reconstruct jets: the di-jet asymmetry [ |133[|134J . Namely, a significant fraction of the di-jet events in 



Pb+Pb collisions at y/s^N = 2.76 TeV shows a large transverse energy imbalance between the trigger 



(or 'leading') jet and the away (or 'subleading') jet. (The left panel in Fig. [42] and the right panel in 
Fig. [2] display such asymmetric di-jet events, as measured by CMS.) One should stress that the criterion 
used to define a 'jet' — the value of the product R = A$ x At/ between the spreadings of the hadron 
yield in azimuthal angle and in pseudo-rapidity — is the same for the leading and subleading jets. 
Moreover, a substantial asymmetry between the two jets exists already in p+p collisions, because of the 
bias introduced by the trigger process (see the histograms in Fig. [43]). But the heavy ion collisions show 
a significant increase in this asymmetry, which becomes more pronounced with increasing centrality. A 
quantitative way to characterize this asymmetry is via the transverse energy imbalance, 

Aj = , (66) 

RL1 + P±2 

where p±i (p±2) is the transverse momentum of the leading (subleading) jet. The normalization in 



Eq. ( [66] ) is useful for removing uncertainties in the overall jet energy scale. Fig. [43] (left panel) shows 
the distribution of the Pb+Pb events as a function of Aj, in different bins of centrality: for the most 
peripheral collisions, this is quite similar to the respective distribution for p+p collisions, as shown in 
figure (a). But for the more central collisions, there is an increase in the fraction of events with relatively 
large Aj = 0.3-1-0.4, which significantly exceeds the respective prediction of the PYTHIA Monte-Carlo 
event generator (which neglects the medium effects). This demonstrates that there is additional energy 
loss by the jet, estimated as 20 to 30 GeV, due to its interactions in the medium. But this energy loss does 



not lead to significant angular decorrelations: as visible in the right panel of Fig. [43} the distribution of 
the subleading jet is still peaked at A$i2 = 7r, like in p+p collisions (and in good agreement with the 
PYTHIA simulations). This implies that the additional energy loss is not to be attributed to rare, hard, 
emissions (which would typically lead to 3— jet events). A careful analysis of the background around the 
away side jet allowed one to establish that the missing energy is in fact associated with relatively soft 



particles (pj_ < 2 GeV) emitted at large angles with respect to the jet axis [ |134[ . This rises the question 



about the physical mechanism which is responsible for such in-medium emissions at large angles. 

As we shall now explain, a natural mechanism in that sense exists in QCD at weak coupling: this is 
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medium-induced gluon radiation, that is, the emission of gluons stimulated by the interactions between 
the partons composing the jet and the constituents of the medium. The most interesting situation — 
originally studied by Baier, Dokshitzer, Mueller, Peigne, and Schiff [ |135[|136[ and independently by 



Zakharov ]137||138[ , following pioneering work by Gyulassy and Wang (139) (see also ]140f[146) ) 



is when the medium is so dense that the gluon formation time is much larger than the mean free path 
of a parton propagating through the medium. In that case, there are many collisions which coherently 



contribute to the emission of a single gluon (see Fig. [44]), leading to a suppression of the radiation 
spectrum as compared to the (Bethe-Heitler) spectrum that would be produced via independent multiple 
scattering. This suppression is known as the Landau-Pomeranchuk-Migdal (LPM) effect. 

The 'gluon formation time' is the typical time that it takes in order to emit a gluon with a given 
kinematics. This concept is quite similar to the 'fluctuation lifetime' introduced in Eq. ([5]), but it is 
instructive to present here an alternative derivation for it, which is adapted to the problem at hand. To 
be specific, consider the emission of a gluon with momentum k and energy uo = \k\ by an energetic 
quark propagating along the z axis. Even though the qqg vertex in QCD is local, the emission process 
is truly non-local, as it takes some time for the emitted gluon to lose coherence w.r.t. its parent quark. 
Namely, when the gluon starts being emitted, its wavefunction is still overlapping with that of the quark, 
so the two quanta cannot be distinguished from each other. But with increasing time, the gluon separates 
from the quark and their quantum coherence gets progressively lost. When the quark is very energetic, 
the gluon is typically emitted at a very small angle 9 ~ k±/uo <C 1 and the coherence between the two 
quanta is measured by their overlap in the transverse space. The gluon is considered as being 'formed' 
(or 'fully emitted') when its transverse separation b± ~ 9 At from the quark becomes larger than its 
transverse Compton wavelength = l/k± ~ l/(uo9). This condition is satisfied after a time 

2uo 2 

At form ~ jy ~ —2 • (67) 
kj_ uo9 2 

(The factor of 2 in the numerator is conventional.) The above argument is completely general: it holds 
for gluon emissions in the medium or in the vacuum. What is different, however, in the two cases is the 
mechanism causing the radiation and the associated gluon spectrum. 

To better appreciate this difference, remember first that an on-shell quark cannot radiate: it can 
produce virtual fluctuations and thus develop a partonic substructure, as discussed in Section [3j but 
energy-momentum conservation prevent these quanta to become on-shell, and hence to separate from the 
parent quark. For the radiation to be possible, the quark and/or the emitted gluon must suffer additional 
interactions, which provide the energy deficit. 

Consider first the situation in the vacuum : the quark is produced in an off-shell state via a hard 
scattering and then evacuates its virtuality via bremsstrahlung. Namely, it emits a gluon with energy uo 
and transverse momentum k± within a time interval At ~ 2uojk\ after the original scattering. These 
values uo and k± are arbitrary (subjected to energy-momentum conservation) and independent of each 



other. However the bremsstrahlung spectrum, Eq. (16), favors the emission of soft (x < 1, or relatively 



small uo) and nearly collinear (9 — » 0, or small k±) gluons, for which the formation time is long. 

The situation is very different for a quark propagating through a dense medium : the quark under- 
goes collisions with the medium constituents, with a typical distance £ (the mean free path) between two 
successive collisions. Any such a collision provides a small acceleration, thus allowing the quark to ra- 
diate. Accordingly, the initial virtuality of the quark is not essential anymore: the quark can now radiate 
anywhere within the medium, and not only within a distance ~ Atf orm after the original hard scatter- 
ing. This implies that the phase-space for in-medium emissions is enhanced by a factor L/Atf orm with 
respect to emissions in the vacuum. Here L is the longitudinal extent of the medium as crossed by the 
quark and is typically much larger than Atf orm , as we shall see. Moreover, the emission mechanism and 
the associated formation time are influenced by the gluon interactions in the medium, which destroy the 
coherence between the gluon and the parent quark and thus facilitate the radiation. 
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Fig. 44: Right: A cartoon illustrating medium-induced gluon radiation by an energetic quark propagating through a 
dense QCD medium. Both the quark and the emitted gluon undergo multiple scattering off the medium constituents 
(represented by the black blobs). Left: A quark-antiquark antenna emitting a gluon (here, from the antiquark leg) 
within the medium with size L. The interactions with the medium are not explicitly shown. 



To estimate Atf orm for medium-induced emissions, we also need the rate at which the (virtual) 
gluon accumulates transverse momentum via rescattering in the medium. We shall later check that 
the successive collisions proceed independently from each other and thus provide transverse momenta 
which are randomly oriented and add in quadrature. This implies that the average transverse momentum 
squared grows linearly with time: ~ qAt, where At is the lifetime of the virtual gluon, as measured 
from the emission vertex, and q is a medium-dependent transport coefficient known as the jet quenching 
parameter, to be specified later. Hence, during its formation, the gluon acquires a typical transverse 
momentum squared k 2 ~ qAtf orm via scattering within the medium. On the other hand, the condition for 
quantum decoherence requires the relation ( [67] ) between Atf orm and k 2 . Together, these two conditions 
determine both the formation time and the typical transverse momentum of the gluon at the time of 
emission: 

2uj 



At form ~ -j & k 2 f ~qAt iorm Atfonn-W— & k) ~ (2uq) 1 ' 2 . (68) 




We thus see that, for medium-induced radiation, k± and uj are not independent kinematical variables 
anymore: the transverse momentum of the emitted gluon is acquired via interactions in the medium 
during the formation time which grows with the energy like Atf orm oc uj 1 / 2 . 

This mechanism for gluon production is operational provided the formation time is much larger 
than the mean free path I, but smaller than the size L of the medium which is available for the emission 
process (the distance traveled through the plasma by the parent quark): 

t < At form < L uj m[n = X -qt 2 < uj < cj c = ^qL 2 . (69) 

These arguments imply that the typical emission angle at the time of formation, Of ~ kf/u, cannot be 
arbitrarily small: 



uj \0J 3 J 



/ — — I ~| ) =^ C EE — — < Of « # m ax = -= • (70) 



Unlike bremsstrahlung, the in-medium radiation does not favour collinear radiation. In fact, the smaller 
is the gluon energy u, the larger is the emission angle Of and the shorter the emission time Atf orm . 
The final spectrum favors indeed the emission of relatively soft gluons with uo <C oj c , for which the 
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formation time is much smaller than the size of the medium, Atf orm <C L, and the emission angle is 
quite large: Of » 9 C . Moreover, after being emitted, the gluons keep interacting with the medium and 
thus get deflected at even larger angles: their average transverse momentum can rise up to a final value 
(kj_) ~ q(L — to) ~ qL, where to is the time at the emission vertex. This phenomenon is known as 
transverse momentum broadening. 

The above considerations show that the medium-induced radiation is very efficient in broadening 
the jet energy in the transverse plane, via the emission of soft (u <C oj c ) gluons, in qualitative agreement 
with the LHC data for di-jet asymmetry. On the other hand, the Raa data for 'high-p^ suppression' 
are probably more sensitive to the emission of harder gluons, with uj ~ uj c , which dominate the energy 
loss by the 'leading particle' (the parton which has initiated the jet). Accordingly, the total energy loss 
AE ~ u c ~ qL 2 scales like the square of the medium length L, and not like L (as one would expect for 
a mechanism where the energy is lost locally, say via elastic collisions in the medium). The reason for 
this scaling with L 2 is, of course, the fact that the actual mechanism at work is non-local : it takes a time 
Atf OTm (uj) to emit a gluon and for uj ~ uj c this time is of the order of L. Since moreover the emission 
can be initiated at any point within L, the overall energy loss scales like L 2 . Reversing the argument, 
one concludes that the stopping length for a particle which loses all its energy inside the medium scales 
like L st0 p ~ E 1 / 2 , where E is the initial energy of that particle. 

For more quantitative studies and applications to phenomenology, one still needs an estimate for 
the jet quenching parameter q. Let us assume, for definiteness, that the medium is a quark-gluon plasma 
with (local) temperature T and weak coupling. The typical collisions which matter for the problem at 
hand are soft, in the sense that the typical transferred momentum is of the order of the Debye mass 
vfiB ~ gT introduced in Section [4] This is so since the soft collisions occur much more frequently than 
the hard ones: the corresponding cross-section is not the transport cross-section evaluated in Eq. ( [56] ) 
(that was the cross-section for scattering at large angles), but rather the total cross-section 

f ,n dcr f 1/1 • n a l a 2 f d9 a s 1 1 

a tot = J dfi- oc J 40^9^^ r.^J-r.-^-. (71) 

This is dominated by small-angle scattering — the integral over 9 is cut off at 9 ~ g by the plasma 
effects — and is larger by a factor l/a s than the transport cross-section ( [56] ). The relevant mean free 
path is obtained by inserting the total cross-section in Eq. (55) : i ~ l/(ncr to t) ~ [Ta s ln(l/a s )] -1 . 
As anticipated, this is parametrically larger than the interaction range 1/mrj ~ 1/ 'gT, meaning that 
successive collisions can be treated as independent. Since on the average there is a transfer Ak\ ~ m 2 D 
of transverse momentum squared per collision, we finally conclude that 

* S ^dT " Hl^s), (72) 

for a weakly-coupled QGP. This is merely a parametric estimate, valid to leading logarithmic accuracy, 
and as such it suffers from the same lack-of-accuracy drawback as discussed for other transport co- 



efficients towards the end of Section |4.3[ it cannot be trusted for phenomenological applications. In 
fact, the whole set-up above described is a bit too idealized to correspond to the actual experimental 
situation. To cope with that, more sophisticated, phenomenological models have been proposed which 
treat the geometry of the collision in a more realistic way (finite volume, longitudinal expansion, time 
and point dependent jet quenching parameter) and include also free parameters. Such models are quite 
successful in describing the data, for both the Raa ratio ( [65] ) and the di-jet asymmetry ( [66] ) (see [147 1 
for a recent discussion and more references), but at the expense of using a rather large (average) value for 
q — considerably larger than the corresponding perturbative estimate to leading logarithmic accuracy. 
This situation is sometimes viewed as an argument in favor of the strong coupling scenario, but it might 
simply reflect the inaccuracy of the current perturbative results. 

Note finally that the theory discussed above has addressed the (medium-induced) emission of a 
single gluon, whereas in reality one expects the in-medium evolution of a jet to involve several succes- 



60 



sive emissions — both by the hard parton which has initiated the jet and by its descendants. (Multiple 
emissions become important when the quantity a s (L/ Atf orm ) — which is roughly the probability for 
one gluon emission — becomes of order one.) Phenomenological models generally assume that succes- 
sive emissions proceed independently from each other, but this is still to be demonstrated: a priori, there 
could be interference effects between emissions by different sources (the various partons forming the 
jet). For jet evolution in the vacuum, one knows that such interference effects are indeed important: they 
lead to angular ordering of the subsequent emissions — the successive emission angles are smaller and 
smaller (2j[l48j. For the case of in-medium radiation, there is so far no explicit calculation of two (or 
more) successive emissions, but there are studies of interference effects in the emission by two sources: 
a quark and an antiquark forming a 'colour antenna' (see the right panel of Fig. [44]) ]146[|149[|150[ . In 



particular, the analysis in ]146[ shows that the interference effects are negligible so long as the antenna 



opening angle (the angle 9 q q in Fig. 44) is much larger than the minimal angle 9 C introduced in ( |70] ). As 
previously explained, the typical emission angles obey this condition already at the time of formation 
and they become even larger at later times, due to the momentum broadening by the medium. This sug- 
gests that successive medium-induced gluon emissions can be effectively treated as independent, thus 
justifying a probabilistic approach to in-medium jet evolution, like in Refs. [151f[153) . 



4.5 The AdS/CFT correspondence: insights at strong coupling 

At several points in the previous presentation, we pointed out observables whose values as extracted from 
the heavy-ion data seem difficult to understand if the coupling is weak, but would be more naturally 
accommodated at strong coupling. These observables include the viscosity-over-entropy ratio r]/s, the 
thermalization time r eq , the jet quenching parameter q, and the heavy quark diffusion coefficient D. In all 
these cases, the hypothesis of a strong coupling must be subjected to caution. First, these quantities are 
measured only indirectly, that is, they are extracted from fits to the data based on complex analyses which 
involve theoretical prejudices (notably, on the overall physics scenario), various assumptions which are 
difficult to check, and a considerable amount of model-building (concerning e.g. the initial conditions for 
hydrodynamics, the geometry of the collision, the theoretical description of multi-particle interactions). 
So, it is fair to say that the systematic uncertainties on these observables are still quite large (even though, 
within a given scenario, they might be strongly constrained by the data). Second, the weak-coupling 
results which serve as benchmarks for comparison are generally leading-order results in the perturbative 
expansion. But, as emphasized in Section [43J the standard perturbation theory (i.e. the strict expansion 
in powers of the coupling) is not reliable for the description of transport phenomena, even if the coupling 
is weak. This is so because of the need to resum finite-density effects like the Hard Thermal Loops 
to all orders, and this has not been done so far for dynamical quantities. The situation becomes even 
more complicated for the far-from-equilibrium situations, as relevant for the phenomenology, where the 
medium effects are not well understood. 

This being said, the hypothesis of a strong coupling is both interesting and intriguing, and not 
easy to refute on the basis of asymptotic freedom or of the current lattice data. Indeed, as already noted 



in Section |4.3[ the QCD coupling g is quite large when evaluated for temperatures a few times T c (the 
critical temperature for deconfinement): g — 1.5-1-2. As also mentioned there, the perturbative series 
at finite temperature is truly an expansion in powers of g (and not of a s = g 2 /47r), so for that purpose 
the coupling is moderately strong. Reorganizing the perturbation theory via appropriate resummations 
of medium effects is one of the possible strategies to cope with this problem. But performing fully non- 
perturbative calculations, whenever possible, is clearly interesting. For thermodynamics, lattice QCD 
is the obvious and pertinent non-perturbative tool. As discussed in Section [43J its results are roughly 
consistent with those of HTL resummations at weak coupling. Yet, as we shall later argue, the lattice 
QCD results for the pressure do not totally exclude a strong coupling scenario. For real-time quantities 
and the non-equilibrium evolution, lattice methods become unapplicable (or, at least, inefficient), so 
it has become common practice to rely on the AdS/CFT correspondence for guidance as to general 
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properties of strongly coupled field theories at finite temperature. (See [154[ for a general review on 



AdS/CFT and Refs. [ |120[|155f[157[ for recent reviews of its applications to a finite-temperature plasma.) 



The AdS/CFT correspondence (or 'gauge/string duality') does not apply to QCD, but to a 'cousin' 
of it, the N — 4 supersymmetric Yang-Mills (SYM) theory, which has a non-Abelian gauge symmetry 
with the 'colour' group SU(iV c ), like QCD, but also additional global symmetries (notably, supersym- 
metry), which strongly constrain the dynamics. These additional symmetries ensure that the conformal 
invariance of the classical Lagrangian is preserved after including quantum corrections — meaning that, 
unlike in QCD, the coupling is fixed and there is no confinement. Accordingly, this theory has probably 
little to say about the zero-temperature, hadronic, phase of QCD, where the non-perturbative aspects of 
QCD are controlled by confinement. Moreover, this is probably not a good model for the QCD dynamics 
in the vicinity of the deconfinement phase transition, where the running coupling effects are known to 
be important, as shown by the lattice results for the 'trace anomaly' in Fig. [37] (left). However, as also 
manifest in that figure, the relative 'anomaly' (e — SP)/e decreases very fast with increasing T above 
T c and becomes unimportant (smaller than 10%) for T > 2T C ~ 400 MeV. Hence, there is a hope that, 
within the intermediate range of temperatures at 2T C < T < 5T C , which is the relevant range for heavy 
ion collisions at RHIC and the LHC, the dynamics in QCD may be at least qualitatively understood by 
analogy with the Af = 4 SYM theory at strong coupling. 

Specifically, the AdS/CFT correspondence is a duality, that is, an equivalence between two theo- 
ries which a priori look very different from each other: (i) the Af — 4 SYM gauge theory mentioned 
above (the 'conformal field theory', or CFT) and (ii) a special, 'type II B', string theory, leaving in 
a curved 10 dimensional space-time with Anti-de- Sitter (AdS) geometr^ This duality is interesting 
in that it maps the strong coupling sector of Af = 4 SYM onto the weak coupling sector of the string 
theory. Accordingly, it allows one to compute observables in the CFT at strong coupling via perturbative 
calculations in the string theory. More precisely, the 'strong coupling limit' to which refers the duality is 
the special limit (g denotes the gauge coupling in Af = 4 SYM) 

A = g 2 N c oo with g 2 < 1 , (73) 

that is, the limit of a large number of colours (N c —¥ oo) taken for a fixed, and relatively small, value of 
the gauge coupling g. This defines indeed a regime of strong coupling (despite g being small) because 
when N c » 1 the effective coupling in the gauge theory is the 't Hooft coupling A = g 2 N c . (This is also 
true for QCD with colour group SU(iV c ).) For instance, the perturbation theory in the multi-colour limit 
is dominated by planar Feynman graphs which are such that each additional loop brings a factor of A. As 
long as A <C 1, the large- N c limit of the theory can be studied in a perturbative expansion in powers of 
A. In the opposite limit A — >> oo (but with g <C 1), one can rely on the AdS/CFT correspondence. In that 
limit, the dual string theory reduces to 'supergravity' (or SUGRA) — a classical field theory in a curved 
space-time with 10 dimensions. From the solutions to the classical equations of motion (e.g., Einstein 
equations), one can unambiguously construct, via the AdS/CFT dictionary, the correlations in the Af = 4 
SYM theory at infinitely strong coupling. 

The Af = 4 SYM theory at finite temperature and A > 1 provides a model for the strongly- 
coupled quark-gluon plasma (sQGP). The corresponding string-theory dual is obtained by adding a 
black hole into the AdS space-time. This is somewhat natural, since, as we know from Hawking, a black 
hole has entropy and generates black-body radiation, so in that sense it behaves indeed like a thermal 
system. The entropy of a black hole is proportional to the area of its event horizon. (No information 
can escape from the volume inside the horizon, so this volume cannot contribute to the entropy.) The 



More precisely, this 10 dimensional space-time is the direct product AdS 5 
de-Sitter space-time, with constant negative curvature, and S 5 is the 5-dimensional sphere, with constant positive curvature; 
see e.g. |154| for details. 
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corresponding, Bekenstein-Hawking, formula can be adapted to supergravity, to yield 



_ Horizon area _ S B h _ * 2 2 3 _ 3 

~~ 4Gio = ~ ~2 -4 S0 ' (74) 

where Gio is Newton constant in 10 dimensions and V3 is the volume of the physical 3-dimensional 
space (as usual, we set h = c = ks = 1). The last equality in Eq. ( |74| ) shows that the entropy density 
s of the J\f = 4 SYM plasma at infinitely strong coupling is 3/4 of the corresponding quantity sq at 
zero coupling ! Hence, in spite of the interactions being so strong, the entropy does not deviate strongly 
from that of an ideal gas. The first correction to this result at strong coupling, of order 1/A 3 / 2 , is also 
known [154[ and it is positive — meaning that the NLO result for s is even closer to the respective 
Stefan-Boltzmann limit. Even though such results cannot be directly applied to QCD, they nevertheless 
suggest that the relatively small deviations — about 20% in the temperature range relevant for HIC's, cf. 
Fig. [36] — between the lattice results for the pressure in QCD and the respective ideal gas limit are not 
necessarily in contradiction with a strong coupling scenario. 

We have previously mentioned that one important prediction of AdS/CFT is the limiting value ( [6T] ) 
for the ratio 77/5, which has been conjectured to be a lower bound of nature (9TJ (as it holds at infinitely 
strong coupling for all the gauge theories having a gravity dual). Let us sketch here the derivation of this 
result (90}[9TJ. As explained around Eq. ( |59| ), the shear viscosity 77 describes the response of the plasma 



to 'shear forces' (its ability to transfer momentum p x along the y direction). This is made precise by the 
Kubo formula which expresses 77 as a 2-point function of the shear tensor T xy in thermal equilibrium: 

r] = lim -L [ dtd 3 xe~ iujt ([T xy (t,x),T xy (0,0)]) T . (75) 

This representation is exact, so in particular it holds at strong coupling. In that case, the string-theory 
dual of the 2-point function in the r.h.s. is the cross-section for the absorption of a soft AdS graviton (the 
supergravity field dual to the energy-momentum tensor in the CFT) by the black hole. This cross-section 
is known from general relativity: it is proportional to the area of the event horizon, like the entropy. Thus, 
in this context, 77 and s are naturally proportional with each other. The proportionality coefficient can be 
explicitly computed, with the result that 77/5 = 1 / (47r) . It is remarkable that the heavy ion data seem to 
favour a value which is close to this conjectured lower bound. 



From Section 4. 1 we recall that the shear viscosity enters the equations of hydrodynamics at lin- 



ear order in the gradient expansion. This corresponds to the fact that, in the respective Kubo formula 



( [75] ), 77 is extracted from the term linear in uj in the small frequency expansion of the 2-point function of 
the shear tensor. By going up to the second order in this expansion, one can similarly extract the trans- 
port coefficients for the second-order formalism (which, we recall, are essential to provide a consistent 
formulation of relativistic hydrodynamics). Interestingly, the calculation of these coefficients turns out 
to be simpler at strong coupling, where one can rely on AdS/CFT for that purpose, than at weak cou- 
pling, where the respective calculations (say, using kinetic theory) require the resummation of infinitely 
many Feynman graphs. And as a matter of fact, the general structure of the second-order terms in the 
equations of hydrodynamics has been clarified only recently, via AdS/CFT calculations at strong cou- 
pling ]158||159[ . Thus, thanks to AdS/CFT, one can study the emergence of hydrodynamics from the 



underlying fundamental field theory in a controlled way (at least on the example of J\f = 4 SYM). 

Since the string theory dual of a finite-T plasma is a black hole in AdS 5, it is natural that the 
phenomenon of thermalization at strong coupling can be studied as the emergence of an event horizon in 
the solution to the Einstein equations. We more precisely mean the equations describing the deviation in 
the AdS 5 metric generated by some matter distribution, which is initially out of equilibrium. For instance, 
the collision between two heavy ions can be modeled as the scattering between two gravitational shock 
waves — two Lorentz contracted shells of matter which propagate against each other and scatter via 
gravitational interactions. This problem has been addressed within the supergravity context, via analytic 
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approximations [160f[162| and via exact, numerical, calculations [163] , with the conclusion that the 
system evolves rather fast towards a (locally) isotropic distribution. Interestingly, it appears that the only 
acceptable solution to the Einstein equations which is boost invariant is the one describing the emergence 
of perfect hydrodynamics (in the sense of Eq. ( [58] )) at asymptotically large times |164[|165| . 

Another phenomenon which is interesting to study at strong coupling is jet quenching — the 
energy loss by a 'hard probe' (energetic parton) propagating through a strongly coupled plasma. The 
corresponding AdS/CFT calculations have been performed both for a very heavy quark (which loses 
only a tiny fraction of its total energy) and for a light parton (quark, gluon, or virtual photon), which can 
be totally stopped in the medium ]120[|156[|157[ (and Refs. therein). Here we shall focus on the second 
case — that of a light, but very energetic, parton with original energy E » T. The corresponding 'dual' 
object on the supergravity side can be a semi-classical string falling into AdSs (in the case of a light 
quark), a pair of such strings (to describe a gluon), or a falling wavepacket carrying the photon quantum 
numbers (for a virtual photon). The respective AdS/CFT calculations [ |166| - |169| revealed that, in all such 
cases, the stopping distance over which the light parton loses most of its energy through interactions in 
the medium scales like 

1/£\V3 



^stop 



T \T J 



(76) 



Note the difference w.r.t. the corresponding result at weak coupling, which in Section [44| has been found 
to scale like E 1 / 2 . This reflects the difference between the respective mechanisms for energy loss, that 
we shall now explain |156[|166[|170| . 

From Section |4.4[ we recall that the mechanism at work at weak coupling is medium-induced 
radiation — the emission of gluons stimulated by the interactions between the radiating system (the 
'hard probe' and its partonic descendants) and the individual constituents of the medium. In general, 
several such interactions can contribute to the emission of a single gluon (the LPM effect), but the role 
of the individual interactions is nevertheless well identified: they provide transverse momentum kicks at 
a rate measured by the jet quenching parameter, Eq. f72| ). This is in agreement with the fact that, at weak 
coupling, the plasma is a collection of elementary constituents, or 'quasi-particles' (cf. Section |4~3] ), 
which are pointlike and quasifree. But at strong coupling, we do not expect such a quasi-particle picture 
to hold anymore — rather, the plasma should look homogeneous, without any microscopic substructure. 
And indeed, the AdS/CFT results like Eq. ( |76| ) can be understood by assuming that the plasma acts on 
the external probe with a uniform force Ft ~ T 2 . This is like a gravitational force in the sense that 
it is fully determined by the local energy density ~ T 4 in the plasma, irrespective of its microscopic 
nature. The effect of this force on a virtual parton (the 'hard probe') is to stimulate gluon emission, via 
medium-induced parton branching [156, 1 66} 170 |. 

Specifically, a partonic fluctuation with energy uo and virtuality Q can decay under the action of 
the plasma force Ft provided the mechanical work W = L Ft furnished by this force over a distance 
L of the order of the lifetime of the fluctuation (L ~ uo/Q 2 ) is large enough to compensate the parton 
virtuality. This condition implies 

1 /Cjy/3 

T 



Q 2 



Q => Q = Qs(w) ~ (ojT 2 ) 



2\l/3 



Q2M 



t) 



(77) 



in agreement with Eq. ( |76| ). More precisely, the above argument provides the typical distance L for the 
occurrence of one branching, but this is of the same order of magnitude as the overall stopping distance; 
indeed, the subsequent branchings involve gluons which are softer and softer, and thus proceed faster 
and faster. For a given energy uo » T, any parton with initial virtuality Qo < Q s can decay in this way, 
including the space-like photon exchanged in DIS (cf. Section 3.1). Accordingly, the quantity Q s (uo) 
plays also the role of the saturation momentum for the finite-T plasma at strong coupling. However, 
unlike at weak coupling, where the phenomenon of saturation requires large gluon occupation numbers 
n ~ l/a s , cf. Eq. (23 ), at strong coupling one can argue that it occurs for occupation numbers of order 
one (B6}[17TJ. 
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In the previous discussion we have implicitly assumed the plasma to be infinite (or, at least, much 
larger than the stopping distance (f76]>). This means that, on the supergravity side, one has studied the 
propagation of the 'dual' objects in a metric describing a black hole into AdSs. The corresponding 
calculations for a finite-size medium are more difficult, in particular because the corresponding metric is 
more complicated. But one can at least heuristically revert the logic leading to Eq. ( [76] ) and conclude that, 
if an energetic parton propagates through the medium over a finit e dis tance L without being stopped, then 
the amount of energy lost by the particle scales like AE ~ L 3 [ 170| . Interestingly, this scaling appears 



to be supported by some of the RHIC data [172J . This result at strong coupling sho uld b e contrasted with 
the corresponding scaling-law at weak coupling, namely AE ~ qL 2 (cf. Section 4.4). This difference 



reflects the fact that the medium-induced parton branching is not a local phenomenon (unlike transverse 
momentum broadening at weak coupling), but is delocalized over a distance of the order of the lifetime 
E/Ql of the decaying parton, which in turn is commensurable with its stopping distance. 

The above picture of medium-induced parton branching can also explain the AdS/CFT results 
for the energy loss and the transverse momentum broadening of a heavy quark propagating through a 
strongly coupled plasma [ |120[|157J (and Refs. therein). In that case, the variables uo and Q which appear 



in Eq. ( [77] ) refer to any of the quanta emitted by the heavy quark: among all the virtual fluctuations of 
the latter, the only ones which can decay (and thus take away energy and momentum) are those which, 
for a given energy uo, have a relatively small virtuality Q < Q s {<jo). (Quanta with Q » Q s {u) cannot 
significantly interact with the plasma and hence they are reabsorbed by the heavy quark.) The energy 
loss is dominated by the most energetic among the emitted quanta — those having a boost factor uo/Q 
comparable to that (denoted as 7) of the heavy quark. These two conditions, uo/Q ~ 7 and Q < Q s (cj), 
imply the following upper limits on the energy and transverse momentum that can be taken away by one 
emitted parton: uo < uj max with u max ~ iQ s {u max ) ~ 7 3 / 2 T and, respectively, Ak± < Q s {u max ) ~ 
7 1 / 2 T. These maximal values control the energy loss and the transverse momentum broadening of the 
heavy quark. By also taking into account the typical duration uo/Q 2 of an emission, one finally deduces 
the following expressions 



dE r— uo 
— — — v A 



dt {u/QD 



V\Ql ~ Vx-fT 2 . (78) 



d<*i) VXQ 2 S ^Q 4 S 



~ V\J^T 3 , (79) 



dt (co/Ql) 7 Q 

for the respective rates. The factor yA in the r.h.s.'s of these equations appears because the heavy quark 
is a semi-classical object which acts as a colour source with a strength of order y/X at strong coupling 
— meaning that it emits a number of quanta (with given uo and Q) of order yfX during the formation 
time uo/Q 2 of one such a quanta. As anticipated, Eqs. ([78])— ([79]) agree at parametric accuracy with the 
respective results of the AdS/CFT calculations |120H157| . 
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